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Abstract 

We use an inverse scattering approach to study multi-peakon solutions 
of the Degasperis-Procesi (DP) equation, an integrable PDE similar to the 
Camassa-Holm shallow water equation. The spectral problem associated 
to the DP equation is equivalent under a change of variables to what we 
call the cubic string problem, which is a third order non-selfadjoint gener- 
alization of the well-known equation describing the vibrational modes of 
an inhomogeneous string attached at its ends. We give two proofs that the 
eigenvalues of the cubic string are positive and simple; one using scatter- 
ing properties of DP peakons, and another using the Gantmacher-Krein 
theory of oscillatory kernels. 

For the discrete cubic string (analogous to a string consisting of n 
point masses) we solve explicitly the inverse spectral problem of recon- 
structing the mass distribution from suitable spectral data, and this leads 
to explicit formulas for the general n-peakon solution of the DP equation. 
Central to our study of the inverse problem is a peculiar type of simul- 
taneous rational approximation of the two Weyl functions of the cubic 
string, similar to classical Pade-Hermite approximation but with lower 
order of approximation and an additional symmetry condition instead. 
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The results obtained are intriguing and nontrivial generalizations of clas- 
sical facts from the theory of Stieltjes continued fractions and orthogonal 
polynomials. 

Keywords: Inverse problem, peakons, Weyl function, cubic string, Pade 
approximation. 

MSC2000 classification: 35Q51, 34K29, 37J35, 35Q53, 34B05, 41A21. 
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1 Introduction 

The PDE 

u t - u xxt + (6 + l)uu x = bu x u xx + uu xxx , (x,i)eR 2 , (1.1) 

is known to be completely integrable if and only if b = 2 or b = 3. The case 6 = 2 
is the well-known Camassa-Holm (CH) shallow water equation 0j. The case 
b = 3 is the Degasperis-Procesi (DP) equation, found by Degasperis and Procesi 
\7\ to pass the necessary (but not sufficient) test of asymptotic integrability, and 
recently shown by Degasperis, Holm and Hone [SJ|S] to be integrable indeed. 

All equations in the family (|l.lfl admit (in a weak sense) a type of non- 
smooth solutions called multi-peakons (peakon = peaked soliton). These take 
the form of a train of peak-shaped interacting waves, 

n 

u(x,t) =^m i (t)e-l x - x *WI, (1.2) 
i=i 

where the positions %i (t) and the heights (t) are determined by the following 
system of nonlinear ODEs: 

n 

<=1 „ (1-3) 
TO fe = (6- l)^m fc m l sgn(x fc - Xi) e~\ Xk ~ x ^ , 
i=i 

for k — 1, ... ,n. Here we use the convention that sgnO = 0, and dots denote ^ 
as usual. Throughout the paper, n will be fixed but arbitrary. 

In the integrable CH and DP cases, explicit solution formulas for Xi(t) and 
nii(t) can be found using inverse scattering techniques, which makes it possible 
to analyze the peakon interactions in great detail. This was shown for the CH 
case b = 2 by Beals, Sattinger and Szmigielski 3 2. The DP case 6 = 3 was 
briefly treated in a short note by the present authors |15j . where we outlined 
a recursive procedure to determine solution formulas for any given n. It is our 
purpose here to continue this study in more detail. Among other things, we 
will give the promised closed form n-peakon solution of the DP equation (see 
Corollary TFZty . 
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This paper is divided into two major parts. In Section |2 we study the 
peakon solutions of the DP equation. First we prove certain scattering properties 
of the peakons, using only elementary methods and the governing equations 
themselves. Then we describe the inverse scattering procedure, based on the 
Lax pair given in [S] . The spectral problem associated to the DP equation is of 
third order and non-selfadjoint. To prove that the spectrum is nevertheless real 
and simple, as well as some other related facts that we need, we make extensive 
use of the scattering properties derived earlier. We state the resulting formulas 
for the general n-peakon solution, and extract precise asymptotic properties as 
t —>■ ±00. 

Sections |3] and 0] are devoted to the forward and inverse spectral problems 
for the discrete cubic string, our main tool for proving the explicit n-peakon 
formulas, but also a very interesting problem in its own right from the point 
of view of operator theory. By the forward cubic string problem, we mean the 
following third order spectral problem: for a given function g(y) > 0, determine 
the eigenvalues z such that nontrivial eigenfunctions 4>(y) exist satisfying 

-4>yyy{y) = * 9(y) <f>(y) for y G (-1, 1), 
0(-l) = <^(-l) = 0, 0(1) =0. 

This spectral problem, which is equivalent under a change of variables to the 
one appearing in the DP Lax pair, is a non-selfadjoint generalization of the 
well-known (selfadjoint) string equation 

<f> V y{y) = z g{y) <t>{y) for ye (-1,1), 

0(-i) = o, Hi) = o, 

describing the vibrational modes of a string with nonhomogeneous mass density 
g(y), attached at the ends y = ±1. The discrete case arises when g — 9i $yi 

is a discrete measure ("point masses gi at the positions $/,•"), so that the eigen- 
functions are piecewise linear in y for the ordinary string and piecewise quadratic 
for the cubic string. 

From what we proved about peakon scattering, it follows that the eigenval- 
ues of the discrete cubic string are positive and simple. We also provide an 
independent proof of this, valid not only in the discrete case, using the theory 
of oscillatory kernels developed by Gantmacher and Krein ^2] ■ 

The discrete (ordinary) string plays a crucial role in finding the general 
n-peakon solution for the CH equation j3||2]. The inverse spectral problem con- 
sists of determining the positions y^ and masses given the eigenfrequencies 
and suitable additional information about the eigenfunctions (encoded in the 
spectral measure of the string, or equivalently in its Stieltjes transform, hence- 
forth called the Weyl function). The solution, found by Krein (and outlined in 
Appendix lAlfor the reader's convenience), involves Stieltjes continued fractions, 
a subject with extremely rich connections to various mathematical areas such 
as the classical moment problem, Pade approximation, and orthogonal polyno- 
mials. 
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The corresponding inverse problem for the discrete cubic string, which we 
solve in Section leads to new and very intriguing generalizations of much of 
this. There are two Weyl function, one of which is unexpectedly determined by 
the other (the proof of this "static" fact uses the isospectral deformation of the 
cubic string induced by the DP dynamics). These two Weyl functions admit 
simultaneous approximation by rational functions with a common denominator, 
similar to classical Pade-Hermite approximation of two (independently chosen) 
functions, but with lower order of approximation and instead an additional 
symmetry condition, and also similar to the Pade approximation provided by 
the convergents of the Stieltjes continued fraction in the ordinary string case. 
The study of this new type of approximation problem is the key to the solution 
of the inverse spectral problem. The coefficients in the common denominator of 
the rational approximants are given by ratios of determinants involving certain 
moment-like double integrals of a spectral measure. The evaluation of these 
determinants is significantly more difficult than for the classical case, where the 
well-known Hankel determinants of moments of the spectral measure appear 
instead. 

It may be expected that Krein's solution of the inverse spectral problem in 
the case of a string with a general mass distribution also admits a generalization 
to the cubic string, but we have not endeavored to address that question here. 

2 Peakon solutions of the DP equation 

2.1 The governing equations 

We begin by summarizing some basic facts about the Degasperis-Procesi (DP) 
equation 

u t - u xxt + Auu x = 3u x ii xx + uu xxx . (2.1) 
It is advantageous to write it as a system, 

= m t + m x u + 3mu x , (2.2a) 
m = u — u xx . (2.2b) 

The function m(x,t) can be thought of as a momentum-like quantity, although 
we are not aware of any actual physical interpretation of the DP equation. The 
n-peakon solution 

n 

u(x > t)=J2 m i(t)e- lx - Xi{t)l (2.3) 

i=l 

arises when m is not a function but a finite discrete measure, 

n 

m(x,t) = 2^m t (t)8(x- x t (t)). (2.4) 

i=l 

Then l)2.2bll is satisfied by construction, while H2.2a|) is satisfied in a weak sense 
if and only if the positions x\{t), . . . , x n (t) and the heights mi(i), . . . , m n (t) 
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satisfy the following system, which was given in [5] and is one of our fundamental 
objects of study here: 



x k = rriie 

1 = 1 



(2.5) 

m k = 2 2J m fc mj sgn(a; fc - Xf) e _ ' iKfc_a!< l, 



for fe = 1, . . . , n. 



2.2 Examples of explicit solutions 

Finding the solution of (|2.5() in the case n = 1 is trivial: mi = constant and 
X\ = Xi(0) + m\t. We can write this in a way that anticipates the general 
pattern: 

b 2 

m iW = TT2> a:i(i)=log&i, (2.6) 

where b x {t) = 6i(0)e t/Al , and the constants Ai = l/mi(0) and &i(0) = e Xl(0) 
are determined by initial conditions. 

Also when n = 2, the equations can be integrated by straightforward meth- 
ods if one changes to variables Xi ± xi and mi ± mi as discussed in [5]. In our 
notation, the solution takes the form 

%-^l&2 

xi(f)=log \ +x l , 
AiOi + A 2 02 

x 2 (t) = log(6i + 6a), 

(A161 + A 2 6 2 ) 2 (2.7) 

AiA 2 (Ai6f + \ 2 b 2 + j^hb 2 ) 

(bi+b 2 ) 2 



mi(t) = 
m 2 (t) = 



^l + x^l + j^b^ 



where b k {t) — &/ c (0)e t / Afe , and the constants Ai, X 2 , 6i(0), b 2 (0) are determined 
by initial conditions. 

For 77, > 3 the peakon equations become considerably more involved, and 
so also the explicit solution formulas, if written out in full. However, if one 
introduces suitable notation for certain symmetric functions of b k 's and Afe's, 
the formulas become very succinct; see Definitions 12.181 and 12.191 for notation, 
Corollary 12.231 for the general statement and Example 12 . 241 for the three-peakon 
solution written out in detail. 

Remark 2.1. To write the formulas for x\(t) and x 2 (t) in 12.7fl in the form 
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given in |3|, substitute 



Ai = — , A 2 = — , P = ci + c 2 , 
Cl c 2 

Mo) = =7^_, Ho) ~ W> 



2(ci-ca)' ( Cl _ C2 )Vr' 

There is a small error in : the term 4c 2 P in their equation (6.8) should be 
4ci-P, and vice versa in (6.9). However, their phase shift formulas (6.10) and 
(6.11) are correct, and consistent with the corrected (6.8) and (6.9). 



2.3 A priori estimates 

In this section we derive some useful information about the peakon dynamics 
directly from the governing equations l|2.5() . Although this is instructive in itself, 
our primary reason for doing this is that we will use this information in order 
to show that certain eigenvalues A^ that appear later in the inverse scattering 
analysis are real and simple. 

There is a distinction between peakons (m& > 0) and antipeakons (m^ < 0). 
As indicated by the formulas for n = 1 and n = 2 above, peakons move to the 
right and antipeakons to the left. If both are present, the system (|2.5ll blows up 
after a finite time when the first collision between a peakon and an antipeakon 
occurs (their corresponding m^s diverge to +oo and — oo, respectively). In the 
case of the CH equation, u(x, t) can nevertheless be continued in a natural way 
beyond the time of collision [2]. The behaviour of the DP equation seems to 
be more complicated in this respect. In this paper we will concentrate on the 
pure peakon case (all > 0), except for Section 12.101 where we will make a 
few comments about antipeakons. 

Assumption 2.2. Except where otherwise stated, we assume the following: 

1. Pure peakons: rrik{t) > for k = 1, . . . ,n. 

2. Ordering: xi(t) < x 2 (t) <■■■ < x n (t). 

We call the set where Assumption 12.21 holds the "pure peakon sector" of R 2 " 
and we denote it by V . In other words, 

V = {(x,m) eR 2 " : x x < • • • < x n , all m, > 0}. (2.8) 

For later convenience, we also introduce the definitions 

xq = — oo and x n +i = oo, (2-9) 

with natural interpretations like e x ° — and tanhxo = — 1. 

Proposition 2.3 (Global pure peakon solution). Given initial conditions 
consistent with Assumvtion \2. e A at some time t = to, there exists a unique solu- 
tion to H2.5(l defined for all t € R and satisfying these initial conditions. There 
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are constants C and D such that mfc(i) > C > and Xk+i(t) — Xk(t) > D > 
for all k and all tgR. In particular, the solution remains in the interior ofV 
for all t G R. 

Proof. The right-hand side of (|2.5[1 is clearly uniformly Lipschitz continuous as 
a function of x^s and m^'s on any compact subset of V, and thus the local 
existence and uniqueness follows from the Picard-Lindelof theorem. Moreover, 
since the right-hand side of l|2.5l) is smooth on V , the solution is smooth on the 
domain of its existence. To show global existance, we need to rule out for the 
solution to: (1) intersect in finite time the boundary of V, (2) blow up in finite 
time. 

It can be verified directly from the equations (|2.5|l that Mi = mi + . . . + m n 
and M n — (jj^.=i m kj fllfc=i(l " e Xk ^ Xh + 1 ) 2 j are constants of motion as long 
as the ordering assumption is satisfied, which is at least true for all t in some 
open interval T containing to- (In Theorem 12.101 we show that Mi and M n 
belong to a set of n constants of motion M%, . . . , M n .) 

The pure peakons assumption implies that Mi and M n are positive, and 
that rrikif) < Mi for all k and all t £ T. Using this estimate on all but 
one factor nik in M n , and noting that each factor in the second product is 
< 1 due to the ordering assumption, we get M n < mk(t)M™~ l n_1 . Hence 
?7ifc(t) > M ra /M" _1 > for all ft and all t € T, so the pure peakons assumption 
holds in T. Since by definition of T the ordering assumption holds there, we 
conclude that both the pure peakon and the ordering assumption hold in T. 
Similarly, (l- e *«.-"*+i) 2 > M n /(Mfl"- 2 ) =: c e (0, 1), hence x k+1 (t)-x k (t) > 
— log(l — c) > for each k. Since the distances between peakons are bounded 
from below by a strictly positive constant, the peakons cannot change places. 
Thus for no finite time can peakons reach the boundary of V . 

Finally, no blow-up of of solutions in finite time is admitted because by the 
second equation (|2.5fl . and with the help of < M„/M™ _:L < < Mi, we 
get M n /M^- X < Xj < M v Thus {M n /M^ x )t < x 3 - x 3 (0) < Mit, implying 
that the solution stays bounded at any finite time and the ordering assumption 
holds for all t (that is, T = R). □ 

By direct analysis of l|2.5|l we obtain the following theorem which shows that 
the peakons scatter and asymptotically behave like free particles travelling with 
distinct velocities. This behaviour in mechanical systems, integrable by the Lax 
method, was first observed and analyzed by J. Moser in |16|. 

The results of the next theorem will be superseded by more precise ones in 
Section ESI when we have at our disposal the explicit formulas for the n-peakon 
solution. However, we will use these a priori estimates in order to show that 
certain eigenvalues Ai, . . . , A„ appearing in the solution formulas must be real 
and distinct. 

Theorem 2.4 (Asymptotics). Let {xk(t),mk(t)} be a solution to (|2.5|) sat- 
isfying Assumvtion \2. e A Then the following claims are true: 

1. The limits mfe(±oo) := lim mk(t) exist and are finite. 
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2. Peakons travel to the right and scatter; i.e., ik(t) > 0, and Xk+i{t) — 
Xk(t) - * oo as t — > ±00. 

5. lim ife(t) = mfe(±oo), and ccfe(t) = mfe(±oo) i + O(l) as t — > ±00. 

^. < mi(oo) < 7712(00) < ••• < to„(oo) and < m n {— 00) < ••• < 
rri2(— 00) < mi(- 00). 

Proof. We will only prove the statements concerning t — ► +00, as the case 
f — ► — 00 is completely analogous. To improve readability let dij — \xi — Xj\ 
denote the distance between the ith peakon and the jth. If we impose the re- 
striction that the smaller index always be written first (i < j), we can remove the 
absolute value signs by Assumption 12 . 2| and Proposition 12. 31 dy = Xj — Xi > 0. 
Then the peakon equations (|2.5() read 



2rn fe 



\j<fc / \i>k 

nii e~ 



N (2.10) 

\i<fc / \i>k 



We prove claims 1 and 2 together. To begin with, it follows from (|2.1U|) that 
Xk > since all m& > 0. Now we proceed by induction on k, starting with n 
and descending to 1. The equation for m n is 



^m % e- d ' n > 0, 



2m„ 

2 — 1 

so m n is increasing, and since m„ is bounded from above by the constant of 
motion M\ = TO «> ^ follows that limt^oo m„ exists and is finite. Hence the 
integral in 

m„(i) = m„(0) exp ^ 2 ^ m 4 (s) e - *"W ds^j , 

is convergent. Since the 771$ 's are positive and bounded away from zero, this 
implies that 

poo 

e -d in (s) ds < 1 < i < n - 1. 



The derivative of the integrand e is bounded: \ik\ < m * = by 

l|2~T(^ . so 



dt 



\(x n -Xi)e- din \ < \x n \ + \±i\ < 2Mi. 



This means that e must vanish as t — ► 00, otherwise the integral would be 

divergent. (Suppose u(t) > 0, \u(t)\ < C, and u does noi vanish as t — > 00; say 
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limsup^^ u(t) = e > 0. Then there is an increasing sequence of points tj — > oo 
(which can be taken to lie more than e/4C apart) such that u(tj) > e/2, and 
hence > e/4 in the interval [tj,tj + e/4C] because of the bound on u. This 
forces J Q u(s) ds = 00.) Consequently, the nth peakon scatters away from the 
ones to its left: 

lim di n (t) = 00, i = 1, . . . ,n — 1. (2-11) 

i — >oo 

We have now proved the base step of the induction. When considering the 
fcth peakon, we make the following induction hypothesis about the peakons to 
its right: for all j > k, 

i. m, (oo) = lim mAt) exists and is finite. 

t — >oc 

ii. The jth peakon scatters away from the ones to its left: e~ dij ( s ) < 00 
and lim duit) = 00 for i = 1, . . . ,j — 1. 

t—>oo 

We need to show that the same holds also for j = k. Integration of the equation 
for TOfe in (|2.10() yields 




Since the right-hand side is an increasing function oft, we know that its limit ex- 
ists, even though it might be infinite. By the induction hypothesis the integral on 
the left-hand side has a finite limit. We conclude that mfc(oo) = lim^oo mk(t) 
exists, and this limit is finite because is bounded. Hence the limit of the 
right-hand side is in fact finite, which, using the same arguments as in the base 
step of the induction, implies e _d<fc W ds < 00 and lim^oo dik(t) = 00 for 
i = 1, . . . , A; — 1; that is, the fcth peakon scatters away from the ones to its left. 
This concludes the proof of claims 1 and 2. 

The scattering property just proved implies that all terms in the equa- 
tion for ik in l|2.10|l vanish as t — > 00, except the lone nik- Consequently, by 
l'Hospital's rule, lim^oo Xk(t)/t = lim^oo ik/l — TOfc(oo), which doesn't quite 
prove claim 3 yet, since at this stage we can only draw the following weaker 
conclusion: 

Xk(t) = mfe(oo) t + o(t) as t — + 00. (2-12) 

We will improve the remainder to 0(1) after using (|2.12|) to prove claim 4. 

Proposition 12.31 tells us that mi is bounded away from zero, hence < 
mi (00), and that the ordering x\ < ■ ■ ■ < x n is preserved, hence (in order not 
to contradict l|2.12|l ) 

wi(oo) < 7712(00) < • • • < m„(oo). 
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To prove claim 4, it remains to show that these inequalities are actually strict. 
Suppose, to derive a contradiction, that not all mfc(oo) are distinct. Then there 
is a sequence of two or more consecutive m-fc(oo) that are equal, say 

m _i(oo) < m (oo) = ••• = m(,(oo) < m,b+i(co) (2-13) 

for some a < b. (If a = 1, then m a _i(oo) is not present here, and similarly 
mt,+i(oo) is not present if b = - n.) We will show that this implies the existence 
of a to such that rh a {t) < and m&(i) > for all t > to- Then we will show that 
one can find t\ > t such that m a (ti) < rrib(ti). Since m a will decrease and raj, 
will increase as t grows, it follows that m a (oo) < m;,(oo), which is the desired 
contradiction. 

The case a = 1 is trivial, since it is clear from (|2.1U|) that mi(t) is decreasing 
for all t. So assume 2 < a < n and rewrite the equation for m a by factoring out 
the exponential in the last positive term (i = a — 1): 



2m Q 



e -d a 




Under the assumption l|2.13|l . equation l|2.12|l implies that d _i )0 = ct + o(t) 
with c = m a (oo) — m a _i(oo) > 0, and d a ^ a +i = o{t) (since m a (oo) = 77i a+ i(oo)). 
Hence the exponent in the term corresponding to i = a + 1 is ct + o(t), which, 
since the factor m a +i in front of the exponential is bounded away from zero, 
means that the second sum tends to infinity as t — > oo. The first sum tends to 
zero by scattering (unless a = 2 in which case it is empty), and m a -\ is bounded. 
Consequently, the right-hand side is negative for all sufficiently large t, so there 
is a t a such that m a (t) is decreasing for t > t a . 

Similarly, m n (t) is increasing so b — n is a trivial case. If 1 < b < n — 1, 
rewrite the equation for uib by factoring out the exponential in the first negative 
term (t = b + 1): 



2 m j, 




b+db 



- nib+i — I ^ m,e 

\i=6+2 



The right-hand side is positive for large t, since the second sum is either empty 
(if b = n — 1) or tends to zero by scattering, rrib+i is bounded, and the term 
corresponding to i = b — 1 in the first sum tends to infinity because its exponent 
is ct + o(t) with c = m,b + i(oo) — m(,(oo) > 0. Hence mb(t) is increasing for 
t > tb, say. Now simply let to — max(t a , tb). Then m a (t) < and m&(i) > for 
all t > to, as was to be shown. 
Finally, observe that 



Xb - x a — m b 
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The integral of this from t = to oo diverges, since d a b —* oo by claim 2. 
However, the contributions from the sums in parentheses are finite by the proof 
of claim 2. This means that 

3 

(rribit) — rn a (t)j dt = +oo, 

o 

which cannot hold if the integrand is nonpositive for all t>to- This shows that 
there is a t\ > to such that TOb(ii) — m a (ti) > 0, which completes the proof of 
claim 4. 

Now it only remains to finish the proof of claim 3 by sharpening the es- 
timate o(t) in (|2.12|) to O(l). It will be convenient here to put the absolute 
value signs back in dij = \ii — Xj\ and remove the restriction i < j. Let 
ci = 5 minj^j |mj(oo) — mj(oo)|. By claim 4 we know that c\ > and then 
H2.12fl shows that there is a C2 > such that on the interval t > 0, say, the 
uniform estimate 

mt(t) e- d «W < c 2 e' Clt (i ? j) (2.14) 
holds. From the equations of motion p.lOjl we get 



Xk{t) - m k (oo)t = Xfc(0) + / (xk(s) - mfc(oo)) rfs = 

(mfc(s) - mfe(oo)) cis. 



x fc (0)+ ^jr*m,(»)e-**Wdflj +^( 



The sum in brackets is bounded as t — > oo, since the integrals are convergent. 
(We saw this already earlier in this proof, but with the estimate (|2.14JI we can 
even see that they converge exponentially fast.) So we need only prove that the 
last integral is convergent. The expression for rhk/Znik in (|2.10() contains n — 1 
terms, all of the form (|2.14(l . so we obtain 



\m k (s) - m k (oo)\ < J |mjfe(£)|df< 

< / 2M 1 (n-l)c 2 e- Cl «d£ = 2M 1 (n-l)— e ~ ClS , 

Js Cl 

and the integral of the right-hand side from s = to oo is convergent. Thus 
Jo { m k(s) — mfc(oo)) ds is absolutely convergent, and the proof is finished. □ 

2.4 Lax pair 

The complete integrability of the DP equation manifests itself in the existence 
of an associated spectral problem, as is usually the case in soliton theory. It 
was shown in [2j that the DP equation is the compatibility condition for the 
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following linear system for the wavefunction ip{x, t;z): 

(d x - dl)i> = zrm/), (2.15a) 
ipt=[z- l {l-dl)+u x -ud x ]i). (2.15b) 

(The term z _1 (l — d x ) could be replaced by z (c — d x ) where c is an arbitrary 
constant, but for our purposes c = 1 is the appropriate choice.) 

Indeed, with L = d x — d x + zm and A = z~ x {l — d x ) + u x — ud x , the 
generalized Lax equation 

L t = [A,L]-3u x L (2.16) 

is equivalent to 

m t + m x u + 3mu x — 0, m x = {u — u xx ) x . (2-17) 

Clearly (|2.17|) is satisfied by solutions of the DP equation (|2.2|) . Conversely, any 
solution of (|2.17|l such that m and u — u xx vanish as |a;| — > oo (as is certainly 
the case for peakons) also satisfies the DP equation. 

Out of all possible solutions to H2.15[) . we distinguish the following one in 
particular: 

Definition 2.5 (DP wavefunction). To a given solution u(x, t), m(x,t) of 
the DP equation l|2.2|l . and to a given complex number z, we associate the corre- 
sponding wavefunction ip( x i t; z) satisfying l|2.15|l and the additional asymptotic 
condition 

tfj(x, t; z) ~ e x , as x — + — oo. (2.18) 

(It is not hard to see that 12.1811 is consistent with H2.15fl if u and its derivatives 
decay sufficiently fast as |ar| — > oo; we show the details for the peakon case in 
Section [2 . 61 below. 1 ) 



2.5 The peakon wavefunction ifj(x,t;z) 

For peakon solutions, when m = 2^™ rrii5 Xi is a discrete measure, the wave- 
function t/j of Definition 12.51 can be written down quite explicitly. Since in this 
case the right-hand side of (|2.15a|l equals zero away from the points Xi, the 
wavefunction is a solution of (d x — d x )ip = in each interval (xu,Xk+i)- In 
other words, tp is piecewise given by the expression 

i>(x,t;z) = A k (t;z)e x + B k {t;z) + C k (t; z) e~ x , 
x e (x k (t),x k+ i{t)) (fc = 0, ...,n). 

(Recall our convention that xq — — oo and x n+ \ = +oo.) 

In order to see how these pieces fit together, consider for a moment t and z 
to be fixed and simplify the notation to 

il>{x) = A k e x + B k + C k e- X , x e (x k , x k+1 ). (2.20) 
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According to l|2.15a[) . ijj xxx is proportional to a Dirac delta at each point x k , 
which means that tp and tp x are continuous while ?p xx has jump discontinuities: 



ip xx (xk+) - ipxx(xk-) = -2zmk4>{xk) (k = l,...,n). (2.21) 
A straightforward computation translates these continuity and jump conditions 
into the following relation between the coefficients in adjacent intervals: 

B k = S k (z) B fc _! , (2.22) 
,Cfe/ \Cfe-i/ 



where the 3x3 matrix S k (depending on z, x k , m k ) is given by 



S k {z) = I - zm k 



1, e" 



(2.23) 



s k (- 



and 



Here / denotes the 3x3 identity matrix. Note that [S'k(z)] 1 
det S k (z) = 1. 

Equation H2.22|) determines tp completely once Aq, Bq, and Co are known. 
The asymptotic condition (|2.18|) is in the peakon case actually an equality, 
ip(x) = e x for x < xi, which corresponds to Aq = 1 and Bo = Co = 0. So for 
peakons we get the following special case of Definition 12.51 

Definition 2.6 (DP peakon wavefunction). To a given solution {x k (t) , m k (t)}^, =1 
of the DP peakon equations l|2. 5J) . and to a given complex number z, we associate 
the wavefunction ip(x, t; z) given by (|2.19l) with coefficients 



f A k (t-z) 

B k {t;z) ) =S k (z)S k . 



f A a {t;z) 
B (t;z) | = 
,C (t;z) 



i(z) ■ ■ ■ S 2 {z) Si(z) 




(2.24) 



(fc = l,...,n), (2.25) 



where the dependence on £ comes from the dependence of the matrices S k on 
x fe (£) and m k (t). 

It is clear from l|2.25|l and the definition of S k that A k , B k , C k are polynomi- 
als in z of degree k, with coefficients depending on all Xi's and m^'s with i < k, 
and with constant terms A^(i; 0) = 1, B k (t; 0) = C k (t; 0) = (since 5fc(0) = I). 
In fact, it is not hard to see that 




k 

E 

p=i 



p-i 



E (iH(na 



1 

-2e x ' 



(-^) P , (2.26) 
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where i} 1 '^) is the set of all p-element subsets I = {i% < • • • < i p } of the integer 

interval [1, k] — {1, .... k}. The empty product (JKZj. m the case p = 1) is 
interpreted as having the value 1. 

Example 2.7. From 12.26fl one obtains 

Ai = 1 — miz, Bi = 2mie Xl z 7 C\ = —mie 2xi z, 

A 2 = l- [m 1 +m 2 ]z+ [mim 2 (l - e* 1 ^ 2 ) 2 ] z 2 , 

B 2 = 2[ mi e Xl + m 2 e X2 ]z - 2[ mi m 2 {l - e Xl ~ X2 ) 2 e X2 ] z 2 , 

C 2 = -[m ie 2xi +m 2 e 2x2 ]z+ [ mi m 2 (l - e Xl ~ X2 ) 2 e 2x2 ]z 2 , 

and 

A3 = 1 — [mi +m 2 + 7713] z 

+ [mim 2 (l - e Xl - X2 f + mim 3 (l - e 21 "^) 2 + m 2 m 3 (l - e X2 ~ X3 ) 2 ]z 2 
- [miro 2 m 3 (l - e :El - 2:2 ) 2 (l - e^" 213 ) 2 ] z 3 . 

Remark 2.8. Although A 2 , for instance, has the same form for all n consid- 
ered as a function of {xi,x 2l mi 1 m 2 }, the coefficient A 2 {t;z) appearing in the 
rightmost interval x 2 < x < 00 in the two-peakon wavefunction is not the same 
as the coefficient A 2 (t; z) appearing in the interval x 2 < x < X3 in the three- 
peakon wavefunction, since the functions {xk(t), mk(t)} are not the same in the 
two cases. In our presentation we always consider n to be fixed (but arbitrary), 
so hopefully no confusion will arise from this slight notational ambiguity. 

The rightmost interval x > x n will be of particular interest, so we define 

A(t;z)\ (A n {t-z)\ 

B{t-z) = B n {t-z) . (2.27) 
C(t;z)J \C n (t;z)J 

2.6 Time dependence of ip 

Now we take a closer look at the time dependence of the coefficients Ak(t;z), 
Bk{t; z), Ck(t; z) appearing in the peakon wavefunction ip of Definition 12. 61 
We know that ip(x,t; z) must satisfy equation (|2.15bJ) : 

tpt = [z _1 (l - d 2 ) +u x - ud x ] ip, 

where u(x, t) = ^™ rrii(t)e~^ x ~ Xi ^ . We claimed before that our choice 

(A (t;z),B (t;z),Co(t;z)) = (1,0,0) 

is consistent with (|2.15b(l . This is easy to see; in the leftmost interval x < X\(t) 
we have u(x,t) = J 2i m l (t )e x ~ Xl{ - t) (= u x ) and tp(x,t;z) = e x (= ip x = ip xx ), 
hence both sides of (j2.15b|) vanish for x < x\(t). 
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Consider now the rightmost interval x > x n (t), where we have u(x,t) — 
YJl m,{t) e x ^- ' x (= -u x ) and V = A(t; z) e x + B(t; z) + C(t; z) e~ x . Inserting 
this into (|2. 15b|> yields 



A t e x + B t + C t e~ x = 



l-dl 



£mie x «) (l + d x ) 



{Ae x + B + Ce~ x ) 



which proves the following: 

Proposition 2.9. When {xk,Tnk\k=i evolve according to the DP peakon equa- 
tions <|2.5[) . the coefficients in the corresponding peakon wave function ip(x,t:z) 
evolve according to the equations 

A t = 0, B t = --2AM+, C t = -BM+, (2.28) 

z 

■where M+ = X)"=i irne Xi ■ 

In particular, the polynomial A{z) is independent of t, so that its coefficients 
are constants of motion for 1)2. 5[l , If we write 

A(z) = 1 - Miz + M 2 z + (-l)"M„z", (2.29) 

then equation <|2.26|l gives the explicit formula (|2.3U|) for M p below. It is also 
useful to write A(z) — n^ =1 (l — j-), where {A^} by definition are the zeros 
of A(z) (see Section f2.7|) . from which we see that M p is the pth elementary 
symmetric polynomial in {1/Afc}. 

Theorem 2.10. The DP n-peakon equations l|2.5|l admit n functionally inde- 
pendent constants of motion Mi, . . . ,M n , given explicitly by 

m p= e (n^)(n (i " e " 3 "" j+i)2 )' (2 - 3o) 

Ie( [i,»]) iei j=i 

where is the set of all p-element subsets I — {i\ < ■ ■ ■ < i p } of {1, ... , n). 

Example 2.11. The constants of motion in the case n = 3 are 
Mi = mi + m 2 + m 3 , 

M 2 = mi m 2 (l - e Xl -' X2 ) 2 + toito 3 (1 - e Xl - X3 ) 2 + m 2 m 3 (l - e* 2 ^ 3 ) 2 , 
M 3 = TOim 2 m 3 (l - e Xl -' X2 ) 2 {\ ~ e 2 ^ 3 ) 2 . 

(Cf. the formula for A3 in Example 12. 71 ) 
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2.7 The peakon spectral problem 



In Definition 12 . 51 for the wavefunction ip, we viewed (I2.15a|) (for fixed t and z) 
as an initial value problem, with initial data given by the asymptotic condition 
at x = — oo. By imposing constraints both at x = — oo and at x = +oo, we turn 
l|2.15a|) (for fixed t) into a spectral problem instead: given m(x) > 0, find the 
eigenvalues z such that 



1px{x) ~ ^xxx(x) = 


z m(x) 




e x ip{x) — i 


• as a; — 


> — oo 


ipx(x) + ip(x) -S 


• as a: — 


> — 00 


e~ x tp(x) — ! 


• as x — 


> +oo 



(2.31) 



has nontrivial solutions ip( x )- These particular boundary conditions, which 
might seem slightly mysterious at this point, will turn out very natural in the 
setting of the cubic string; see Theorem 13.11 In any case, they take a much 
simpler form in the peakon case, which is what we are interested in here. 

Theorem 2.12. Let m = 2J^™ S Xi be a discrete measure with nrii > and 
x\ < • • • < x n . In this case, the spectral problem (|2.31|l has n distinct positive 
eigenvalues 

< Ai < ■ • • < A„. 

Proof. We saw in the previous section that the solution to (|2.15a|l in the discrete 
case is given by (|2.20|l and (|2.22() . and then the boundary conditions in l|2.31|l 
reduce to 

C =0, B = 0, A n - 0. 

It is clear that we may normalize ip by taking Aq = 1, and then A n (z) is 
exactly the polynomial we considered before; its coefficients depend on the Xk's 
and ink's, and are given explicitly by Ij2.26|l with k = n. Consequently, the 
eigenvalues of Ij2.31|l are in the discrete case simply the zeros of this nth degree 
polynomial A(z) = A n (z). 

Now let {xk,mk} evolve according to the peakon equations l|2.5|l . Even 
though this changes the measure m, the eigenvalues of (|2.31() remain unchanged, 
since the coefficients of A(z) are constants of motion according to Section VTl^ 

As t — + oo, the peakons scatter as described by Theorem 12 .41 the m^'s tend 
to distinct positive values mfc(oo), and xu+i — Xk — > oo. By Theorem 12.101 
this reduces the coefficients of A(z) to the elementary symmetric polynomials 
in {mfe(oo)}, so that A(z) = rifc=i(l — zm k{oo)). It follows that the eigenvalues 
are {l/mfc(oo)}, hence positive and distinct. □ 

Remark 2.13. Since mi(oo) < • • • < m n (oo) and m\{— oo) > • • • > m„(— oo) 
by Theorem 12.41 and since the argument of the above proof holds equally well 
as t —t — oo, it follows that m„ + i_fc(oo) = l/A^ = ?«&(— oo). In other words, 
the same asymptotic velocities occur as t — > +oo and t — > — oo, but in reverse 
order. We will study asymptotic peakon properties of this kind in more detail 
in Section l2. 91 after giving the explicit formulas for the n-peakon solution. 
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It is clear from the isospectral deformation argument in the proof above 
that the numbers {xk,mk} r k l =1 are not uniquely determined by the eigenvalues 
Afe alone. The eigenvalues must be supplemented by further spectral data in 
order for the inverse spectral problem of reconstructing m(x) to have a unique 
solution. 

Definition 2.14 (Extended spectral data {Afc,6fe,Cfe}). Given a point in 
the peakon sector V, 

x\ < ■ ■ ■ < x n and mi , . . . , m n > 0, 

we assign to it a 3n-tuple of numbers {Afe, bk, Ck} k =i y i a the following construc- 
tion. Let the polynomials (A(z), B(z),C(z)) — (A n (z), B n (z),C n (z)) be given 
by jmni with k — n. Let 

< Ai < ■ • • < A„ (2.32) 

be the zeros of A[z) (which are positive and distinct by Theorem 12. 12J) . and let 
bk and Ck be the residues in the partial fraction decompositions 

-B(z) b k C{z) - B{z) ^ c k 

2zA(z) j^z-Xk' 2zA(z) j^*-V 

For future convenience (see Theorem 13 .5|) . also let 

A o = 0, 6 = 1, c = i. (2.34) 

Theorem 2.15. When {xk,mk}^ =1 evolve according to the DP peakon equa- 
tions (|2.5|l . the extended spectral data {Afe, 6fc, Cfc}^ =1 evolve according to the 
equations 

A fe = 0, 6* = -^., c fc =6 fc + V^. (2.35) 
Proof. We have already seen that the eigenvalues Afe are constant. Let 



B{ Z ) " bk 



2zA(z) f~[z-Xk 



and 



c( z ) = c{z) = v Ck 



2zA(z) z- Afe 



k=l 

Equation it^p) shows that A(0) = 1 and B(z) = 2M + z + 0(z 2 ), so that 
w(0) = -M+, where M + = m, fc e x *. The time evolution of A(z), B(z), C(z), 
given by (|2.28|) . now translates into 



B f (z) = B(z)/z - 2A(z)M + co(z) - u(0) 
Wt ^ _ 2zA(z) ~ 2zA(z) - z 



17 



and 

CtO) -B(z)M+ , s , s 

Comparing residues at z — X k on the right- and left-hand sides of these equations 
one obtains at once l|2.35|) . □ 

Corollary 2.16. The c k 's are determined by the X k 's and the b k 's through the 
formula 

c k = X k b k V - 3 — (k = l,...,n). (2.36) 

(Since Cq = 1/2 and bo = 1 6?/ definition, the formula holds trivially also in the 
case k = ; if interpreted as being a limit as Xq — > 0+.) 

Proof. From the second equation in (|2.35(l it is immediate that b k (t) — b k {Q)e l ^ Xk . 
Inserting this into the third equation in l|2.35|) we get 

d, ... , ^ b k (O)b0) ( t t\ 

-(c k (t) b k {t)) = ^ exp ^- + - j . 

From H2.26J) and Theorem 12.41 it is clear that C{z) — > as t — > — oo, hence £(2) 
and its residues — also vanish as t — > — oo. Thus, integrating the ODE 
above from — oo to i yields 

m , ^ 6 fc (0)6 3 -(0) f * t\ ^ A fc b fc (*)*>j(*) 



Using the definitions Ao = and bo — 1, the 6fc on the left-hand side can be 
incorporated in the sum by summing from j > = instead of j = 1, and this gives 

eu. □ 

Theorem 2.17. TTie 's and c k 's are positive. 

Proof. First we show that b k ^ 0. By fe fe = -5(A fc )/2A fc A'(A fc ), so it 

is enough to show that B(X k ) ^ 0. Let ip k (x) denote the eigenfunction cor- 
responding to the eigenvalue X k . Since m has compact support, ip k and all 
its derivatives vanish as x — > — oo, while as x — > +oo the derivatives vanish 
and ^fc(oo) = B(X k ). Thus it suffices to show that ip k (oo) ^ 0. Since ij} k 
satisfies ip x — tpxxx = X k mip we readily obtain, after the standard steps of mul- 
tiplying by integrating and using the boundary conditions, that ^^(oo) 2 = 
Afc J^x, m ^P k dx > 0. This proves that b k ^ for all k. 

To prove b k > we use induction on the number of peakons n. The case 
n = 1 follows from explicit formulas. We assume the claim to be valid for an 
arbitrary configuration of n— 1 peakons. Then we use the following deformation 
argument. 
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From the explicit formula l|2.26[) for A(z) and B(z), we see that by letting 
X\ tend to — oo (keeping the other variables fixed) we can make the coefficients 
in the rational function 

B{z) 



w(z) = — — = > 



2zA(z) fr[z-\k 
come arbitrarily close to the coefficients in the function 

-(1 - zm 1 )B{z) 



2z(l -zm 1 )A(z)' 



(where B and A are the polynomials computed for the (n — l)-peakon con- 
figuration consisting of all but the first peakon), which has a partial fraction 
decomposition 

— hY^ — (ki := l/mi). 

Z — Kl £ — ' Z — Kk 

(If 1 /mi happens to coincide with one of the eigenvalues K2 , . . . , K n of the (n— 1)- 
peakon configuration, then just change m\ a little before x\ is sent to — oo so that 
all Kfc's are distinct.) The poles and residues {Afc,6fc}" of the rational function 
ui{z) depend continuously on its coefficients, which in turn depend continuously 
on {xk,mk}i ■ Since the 6fc's are nonzero, it follows that they cannot change 
sign during this deformation, which brings the Afc's arbitrarily close to the Kk's, 
and b\, . . . , b n arbitrarily close to 0, tfe, • • • , d n . By the induction hypothesis d^'s 
are all positive. This implies that b%, ■ ■ ■ , b n must have been positive to begin 
with. 

Moreover, we know from the proof of Theorem 12.1 51 that the strictly positive 
function M + = Y^k=i m ke Xk satisfies M + — — w(0) = J2k=i bk/^k- Letting the 
peakons evolve according to the DP equation, b\{t) = &i(0)e*/ Al is the dominant 
term in this sum as t — > oo. This leads to a contradiction if b\ < 0, hence b\ > 0. 

Finally, since all Afc's and bk's are positive, equation l|2.36|l shows that all 
Cfc's are positive as well. □ 

2.8 The explicit n-peakon solution 

In the extended spectral data {A^, bk, Cfe})? =1 investigated in the previous section, 
the quantities Ck play only an auxiliary role (although important), since they are 
determined by the Afc's and frfc's. The primary objects are {Afc,6fc}J? =1 , which 
we will refer to as spectral variables. Definition 12.141 gives the prescription for 
going from peakon variables {xk,rrik\ to spectral variables {Afc,£>fc}. We can 
view this as a map 

S :V^TZ 

< \ (\h\ (2 - 37) 
(x,m) i * (A, b), 

where the "pure peakon sector" V is defined by H2.8J1 , while the region of allow- 
able spectral data is 

K = {(A, b) e R 2 " : < Ai < • • • < A„, all b % > 0}. (2.38) 
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Theorems 12.121 and 12.171 show that S does map V into 1Z. To justify the ter- 
minology "spectral variables" , we need to show that S is injective; in fact we 
will show that S is a bijection of V onto 1Z, and hnd 6> -1 explicitly. We will do 
this by transforming the spectral problem l|2.31|l to the equivalent cubic string 
spectral problem l|3.3|l and studying the inverse spectral problem in that set- 
ting. We will state the main result in Theorem 12. 221 in order to explore some of 
its consequences below, but the central part of its proof will have to wait until 
Section IO 

Once we know how to invert the change of variables (|2.37|) , we can also solve 
the DP n-peakon equations l|2.5f) . since in new variables they simply take the 
form of the linear equations i|2.35[) . Using an evolution operator describing the 
solution to (|2.35l) . 

(2 39) 

(A fe A) >- (AfcAe*/**), 
we can write the n-peakon solution as 

(x(t),m(t)) =5" 1 o$ i o5( J :(0),m(0)). (2.40) 

However, finding the formulas for <S -1 is quite a difficult problem, which will 
occupy us for most of the remainder of the paper. We start by introducing 
necessary notation. 

Definition 2.18 (Various notation). For k > 2, let 

A(xi, . . . ,x k ) = Y\_(xt - xj), r(xi, . . . ,x k ) = Y[{x t + Xj). (2.41) 

i<j i<j 

For k > 0, we recall that denotes the set of fc-element subsets I = {i\ < 

■ ■ ■ < i k } of the integer interval [1, n] = {1, . . . , n}. For I <E , let 

A / = A(A ll ,...,A,J, r 7 = r(A il ,...,A i J ) (2.42) 

with the special cases A0 = = A{jj = T^j = 1. For two disjoint ordered 
sets I and J we will also use occasionally 

a Ij= n (^- A i) 2 . r ^,/= n (^+ a j )' ( 2 - 43 ) 

ieijeJ ieijeJ 

with the convention that the symbols equal 1 if one or both sets are empty. 
(Note that Tj u j = FiTjTi t j, and similarly for A 2 .) Furthermore, let 

A/ = ^QAi, &i- = JJbi, 

iei iei 

with the proviso \% = = 1. 

Using this notation we now define the symmetric functions that appear in 
the explicit solution formulas for DP peakons. 
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Definition 2.19 (Symmetric functions). 

A? 



E 



Ie( [i.»]) 



E 

J6( [ij»]) 



A 



7 Ar6 



JO/, 



and 



f/fe+i 



— UkVk — Uk+iVk-i- 



(2.44) 



(2.45) 



(To be explicit, U = V = 1, and f7 fc = Vfe = for fc < or k > n.) 

We can evaluate Wk explicitly in terms of (A, b) as follows. 
Lemma 2.20. 



Wk= E 

Ig( [ V ]) 



E E b M 2 

/nj=0 



,,, f i A / A /,jW 

rvr 



/ i JUJ 



(2.46) 



E 



A^A^r c r D 



V CUD=J 
\C\ = \D\=rn 
min(C)<min(D) 



Proof. For brevity, let = Af/I 1 / and = AJ j/T It j. Consider first the 

contribution from UkVk, which is a sum of terms a&a)(^ B^sbs), with A, £? G 

ft 1 '™]). Since 6a and bs each contain only distinct fej's, we have &a6b = bjbj 

where I — AO B and J = AAB (symmetric difference); note that I fl J = 0. 

The terms for which / = ^4 = _B and J = give rise to the first sum in (|2.46() . 

Otherwise, |/| = fc — m and |J| = 2m for some 1 < m < fc, and the sets E = A\I 

and F = B \ I partition J = E 1 U F into two disjoint sets with \E\ = \F\ = rn. 

With this notation we have (* a&aX* bAb^s) = (*^/ ! ./A/)(* £ ;*_fAf)(6?6j)- 

The contribution from [7fc+iVfc_i is similar, except that A G (fc+i) an d 

[i,«h 
fe-i; 

J = G U H, where G 
= m — 1. Hence, 



£? G (fcl'i) 80 that the case A = B can never happen. In the disjoint partition 

n + 1 and 



A \ / and H = B \ I, we have |G| 



^ = 



fc 

m=i ^(^) 

/nj=0 




E 

GUH=J 
|G|=m+l 
|ff|=m-l 



* G *H A H 
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The expression in brackets equals 1/Tj times the left-hand side of the symmetric 
polynomial identity 



G,H 



/ ^E F Vr / J 

E\JF=J GUH=J 
\E\ = \F\=m \G\=m+l 

\H\=m-l 



2 m+i Aj (A c A D ) 2 r c r fl , (2.47) 



v. CUD=J 
\C\ = \D\=m 
min(C)<min(D) 



substitution of which yields the desired formula 

To prove (|2.47l) we use induction on m = \ J\ /2. Without loss of generality, 
we can take J = {1, 2, . . . , 2m}. When m = 1, both sides reduce to 4 A1A2. For 
m > 1, evaluating l|2.47|l at A 2m _i = A 2m = cyields 4c 2 (n-I™^~ c) 2 (Xj+c)) 
times the left- and right-hand sides, respectively, of the m — 1 case of (I2.47|) . 
which holds by the induction hypothesis. Hence, by symmetry, (|2.47|l holds 
whenever two variables are equal. This means that the difference between the 
left- and the right-hand sides is divisible by Aj, but since that difference is a 
symmetric polynomial it must be divisible by Aj. However, the difference is of 
degree 3m — 1 in Ai, while Aj has degree 2(2m— 1) in Ai. Since 3m— 1 < Am— 2 
if m > 1, we conclude that the difference is identically zero, and (|2.47() is 
proved. □ 

Remark 2.21. Similar polynomial identities will be used later in the proof of 
Lemma T4. 101 see also Appendix IE! The constraint min(C) < min(D) eliminates 
redundancy; if it is removed, then every term appears twice in the sum, and 
2 m+1 should be changed to 2 m . 

Theorem 2.22 (Inverse of S). The change of variables S : V — > 1Z is a 
bijection. Its inverse 

S- 1 :K->V , 



(A, b) i-> (x, m) 



is given by 



1 U k (Uk) 2 (Vfc-1) 2 , . , . 

^fc'=logT7 — ' mk ' = (fc = l,...,n), (2.49) 

Vk-l vVkWk-i 

where k' = n + 1 — k. 

Proof. For the proof, denote the map (|2.49|l by T instead of S^ 1 . In Sect ion l4~4l 
we will prove that T oS — id-p, using the solution of the inverse spectral problem 
for the discrete cubic string. In words: provided that (A, b) = S(x,m) is in the 
range of S, we know that (x, m) = T(A, b). To complete the proof, we need to 
show that the range of S is all of 1Z (that is, all allowable spectral data really 
correspond to a peakon configuration). 
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First note that T maps TZ into the pure peakon sector V . Indeed, from the 
definitions of Uk and 14, and from the explicit formula l|2.46[) for Wk, it is clear 
that all Uk, 14, and Wk (for 1 < k < n) are positive when (A, b) e TZ; then 
<|2.49[1 produces a point (x,m) = T(\,b) with all nij > (obviously), and also 
with x\ < ■ ■ ■ < x n since 



Recall the definition of S: given (x, m) G we form the polynomials A(z) 
and B(z), whose coefficients are Laurent polynomials in {e Xk ,mk} k l = i- Write 
A(z; (x, m)) and B(z; (x, m)) to denote this dependence explicitly. The A^'s are 
defined as the zeros of A(z), and the 6fc's as the residues in — B(z)/2z A(z). 

Now we use the algebraic nature of the mappings involved. The coefficients 
in the polynomial A{z\ T(A, 6)) are rational functions of {A^, 6/c}5J =1 which, since 
T o S = id-p, agree with the coefficients in Ilfc=i(l — z /^k) on the range of S 
(which is an open set in 1Z since S is a homeomorphism) . Hence A(z; T(A, 6)) = 
rifc=i(l — z /^k) identically as a rational function of (A, b); in particular this 
relation holds for all points (A, b) £ TZ. Similarly, the coefficients of B(z; T(A, b)) 
agree with — 2z{%2™_ 1 bi/(z — Ai)) Ilfc=i(l — z /^k) on the range of S, hence on 
all of TZ since they are rational functions. This proves that S o T = id-ji- 
Consequently, the range of iS is TZ, and T = S^ 1 , as claimed. □ 

Corollary 2.23 (The n-peakon solution). The solution {xk(t),mk(t)}^ =1 
of the Degasperis-Procesi n-peakon equations (|2.5|l is given by (|2.49|l with 



where the constants {Xi, &i(0)}™ =1 are determined from the initial conditions 



Example 2.24 (The three-peakon solution). For n — 1 and n = 2, equation 
H2.49|l reduces to (|2.6|l and (j2.7|l . respectively. The solution for n = 3, when 
simplified by taking into account that some of the t/fe's and T4-'s equal or 1, 
takes the form 




h{t) = fo 4 (0)e*/ A * 



(i = l,...,n), 



{x k (0),m k (0)} 



n 



x 3 (t) 



logC/i, x 2 {t) 




Xl{t) 




(2.50) 



m 3 (t) 



(Ui) 2 



m 2 (t) 



(^) 2 (vo 2 



mi(t) 



U 3 (V2) 2 
V 3 W 2 ' 



Wi 



W 2 Wi 
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where, with bi = 6^(0) 



Ui =6i + 6 2 + 6 3 , 

Vi = Xih + A262 + A363, 

Al + A2 Al + A3 A2 + A3 

t r ( A l ~ A 2) 2 . . , , (Al - A3) 2 (A 2 - A3) 2 

V2 = — — — — A1A26162 H — — - — A1A36163 H — — — A2A36263, 

Ai + A2 Ai + A3 A2 + A3 

TJ _ (Ai-A2) 2 (Ai-A 3 ) 2 (A 2 -A 3 ) 2 

3 (Ai+A 2 )(Ai+A3)(A 2 + A 3 ) 1 2 3 ' 

v ( A i- A 2) 2 (Ai-A 3 ) 2 (A 2 -A 3 ) 2 , 

^ = (Ai+A 2 )(Ai+A3)(A2 + A 3 ) AlA2A 3 fo i &2 ^ 

and consequently 

Wi = U-lVi - U 2 

W 2 , , W 2 , 4A 1 A 2 l l 4AlA 3 4A 2 A 3 

= A161 + A 2 6 2 + A363 + - — ——6162 + - — ——6163 + - — —--6263, 

Al + A2 Ai + A3 A2 + A3 

W 2 = U 2 V 2 - U3V1 

(Ai - A 2 ) 4 . . 2 (Ai - A3) 4 2 (A 2 - A3) 4 

= TTTTw AiA 2 (6i6 2 ) + . „ AiA 3 6163 + jt— — , 2 A 2 A 3 (6 2 63 

(Ai+A 2 )^ (Ai + A 3 )^ (A 2 + A 3 )^ 

4A1A2A3616263 / 2 x 

+ (Ai + A2)(Ai+A 3 )(A2 + A 3 ) X (, (Al " A2) (Al " A3) 61 
+ (A 2 - Ai) 2 (A 2 - A 3 ) 2 6 2 + (A 3 - Ai) 2 (A 3 - A 2 ) 2 6 3 

We take the opportunity here to correct a mistake in our previous paper |15| . 
where the denominator W\ in m 3 (i) was incorrectly stated with (Aj + A 3 ) 2 
instead of A; + Xj . 

2.9 Peakon asymptotics and scattering 

We will now extract information from the explicit solution formulas (|2.49|l about 
how the peakons interact, and about their asymptotic behavior as t — > ±00. Re- 
call that we number the eigenvalues so that < Ai < • • • < A n , and consequently 

1 1 

An Al 

Also recall we use k' = n + 1 — k to denote the complementary index to k. 



24 



Theorem 2.25 (Asymptotics). 

x k (t) ~ y- + log6 fc (0) + lQ g A \ i u » as< 

**< (*) ~ f + log 6 fe (0) + £ log tyV^t ' as * ■ 
Afc r-f (A, + Aft )\i 



(2.51) 



-oo, 



and 



wher 



m k (t) = l/\ k + 0(e tSk ), ast^-oo, 
m k ,(t) = l/\ k + 0(e~ tSk ), ast^+oo, 



(2.52) 




k = 1, 

fe = 2,...,n-l, (2.53) 
fe = n. 

Proof. This is simply a matter of identifying the dominant terms among the 
6ft 's. Because of the ordering of the eigenvalues, b\{t) — 6i(0)e*/ Al domi- 
nates as +00 and b n (t) = 6 n (0)e'/ A ™ dominates as t — > — oo. More gen- 
erally, 6162 ■ ■ • 6ft grows exponentially faster as t — > +00 than any other product 
6i 1 6j 2 • • • with ii < ■ • • < ift. It follows that, as < — > +00, 

rr (Ai2...fc) 2 

C/fc - — 6i6 2 • • ■ Oft, 

1 12. ..ft 

and 

Vft ~ (^12. ..ft) ^ & ^ ■ • -6 fe , 

1 12. ..ft 

while (by similar reasoning) for t — > — 00 each index on the right-hand side is 
replaced by its complementary value i' = n + 1 — i. Hence, as t — > +00, 

C4 _ Y A i2-fc V r r2...(fc-i) 6ft = bk ^ e t/\„ tt ( A ifc) 2 



Vfc-i V A i2...(ft-i) / ri2...fc AiA 2 ...Afc_i TjftAi 

Since Uk/Vk-i = exp(x k >), the case t — > +00 of Q2.51.jl follows by taking loga- 
rithms. The case t — > —00 is similar. 

As for Q2.52[) . we have already seen fRemark l2.13|) that m, k (— 00) = 1/Aft = 
77ift'(+oo), and this can of course also be concluded by comparing coefficients 
of the dominant term (61 ■ ■ • 6ft_i) 4 (6fc) 2 (as t — > +00) in the numerator and the 
denominator of the formula for m k > (t). As for the correction term, for 1 < k < n 
a routine computation gives 

1 (i + Aft^(i + £ft^ 

Tfift' (t) = ■; n ; h lower-order terms 

fc (l + Cft^)(l + 73 fe> -M 
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as t — ► +00, where Ak, Bk,Ck, Dk are some i-independent and computable co- 
efficients. This implies 

m k' (t) ~ t — = ctk ~~r~ + Pk 7— — I- lower-order terms 

Afc O fe Ofe_i 

which gives the stated result as i — > +00. A similar analysis applies to the cases 
k = 1 and k — n, where there are no terms involving bk-i or bk+i, respectively. 
Again, the case t — * — 00 is similar. □ 

Corollary 2.26 (Phase shifts). The peakon with asymptotic velocity 1/A& as 
t — » ±00 experiences the phase shift 

lim ( Xk'[t) ~r~ J lim [ Xfc(i) - 

t^oo y Afc/ t-v-oo \ Afe 

fe — 1 /A i\9 71 / \ \\2 



Vlog ( A '~ A fe) _ y W ( A ' Z Afc ) . (2.54) 



Remark 2.27. For every t there is a well-defined notion of what "the peakon at 
site k" means, namely, the peaked wave mfc(i)exp(— \x — Xk{t)\) described by 
nik and Xk- For other types of soliton equations, individual solitons are usually 
only identifiable as they scatter when t — > ±00. In that case, the highest/fastest 
soliton, which is to the far left as t — > —00, reemerges to the far right as t — > +00 
after some interactions during which the identities of the solitons are blurred, 
and similarly for the other solitons. According to this point of view, the "fcth 
fastest/highest peakon" (the one with asymptotic velocity 1/Afe) moves from 
occupying site k as t — > — 00 to occupying site k' as t -> +00, and does not 
really have a well-defined identity during the interactions causing this change. 
It is the phase shift of this fcth fastest peakon that is given in l)2.54[l . 

Remark 2.28. Note that the total phase shift in (|2.54() equals the sum of phase 
shifts resulting from pairwise interactions; since the peakons reverse their order 
(according to the interpretation above), the fcth fastest peakon is overtaken by 
the ones initially to its left (giving the contribution Y] ■ ~. ) and overtakes the 
ones initially to its right (which contributes — X)"=fc+i)- 

We see that the interaction of a faster peakon overtaking a slower one (with 
asymptotic speeds c/ = 1/A/ and c s = 1/A S , respectively, where / < s so that 
Cf > Cg) results in the slower peakon shifting by the amount 

, (A/-A s ) 2 (c/-c s ) 2 
1°§ 7^ = l°g — 



(A/ + A S )A/ (c f +c s )c s 
which is positive iff c/ > 3c s , while the shift of the faster peakon is 

which is always positive. This generalizes the results obtained in |5] for the case 
n = 2. 
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Remark 2.29. The phase shifts obtained in the Camassa-Holm peakon case 
in [2], as well as finite Toda phase shifts obtained by Moser in ^H], have a 
slightly different structure; namely, in the case of CH peakons we have (after 
some rescaling) 

(Aj — Afe) 2 (A; — Afe) 2 



»=1 1 i=k+l ' l 



2.10 A few comments about antipeakons 

So far we have only been discussing the pure peakon case where all m., > 0. 
Everything we have said applies (mutatis mutandis) also to the pure antipeakon 
case where all m« < 0, because of the following natural symmetry of the problem. 

Proposition 2.30 (Antipeakons). Let {xk(t),rrik(t)} be a global solution to 
l|2.5|l satisfying Assumption ^. 'A (pure peakons and ordering). Then 

= -™>n+l-k(-t), x*k(t) = x n+ i- k (-t) (2.55) 

is a global solution to (|2.5I) satisfying m* k (t) < 0, X*(t) > x^t) > ■ ■ ■ > x* n (t). 

Proof. Straightforward. □ 

However, when both peakons and antipeakons are present, the situation 
is more complicated, because of the singularities that occur after finite time 
when a peakon travelling to the right collides with an antipeakon travelling 
to the left. Most of our proofs break down, since they rely on the existence 
of a global solution (Proposition and its scattering properties as t — > ±oo 
(Theorem EH)- 

We suspect that the spectrum is still real and simple in the presence of 
antipeakons, with the same number of positive (negative) eigenvalues as peakons 
(antipeakons). This is easy to prove for n = 1 and 2, and also agrees with the 
analogous situation for the CH peakons |2J. However, in contrast to the CH 
case, the spectral measure /i = bk S\ k need not be positive in the mixed 
peakon-antipeakon case (it is not hard to find counterexamples for n = 2). 

The explicit solution formulas (|2.49(l remain valid as long as they make 
sense, because of their algebraic nature. However, they break down completely 
if some Xi + Aj = 0; for example, in the case of a totally symmetric peakon- 
antipeakon collision (equation (6.13) in |S]). And even if all Ai + Xj ^ 0, the 
solution provided by 12.49fl will in general not satisfy the ordering assumption 
xx < ■ ■ • < x n , or even be defined, for all t. This is because [//., 14, and W% need 
not remain positive (even if positive at t = 0) when some bk < 0. In intervals 
t € {ti,t 2 ) where the ordering assumption is violated, l|2.49|l gives a solution of 
H2.1L)fl (now with some dij = x% — Xj < 0), but this will not be a solution of the 
original peakon equations (|2.5|l containing \dij\ instead of dij. 

All this is also in contrast to the CH case, where the solution formulas for 
Xk(t) and Xk+i{t) always automatically satisfy Xk < Xk+i for all t, with equality 
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only at the instant of a peakon-antipeakon collision where rrik{t) and mk+i(t) 
diverge to ±00. This provides a continuation of the solution past the collision. 
It seems as if continuing DP peakon-antipeakon solutions past collisions would 
require some kind of gluing of different solutions from different time intervals, 
but the details are not clear at present. 



3 The cubic string 

The peakon spectral problem (|2.31(l is equivalent under a change of variables 
to what we have called the cubic string problem, an interesting non-self adjoint 
third order generalization of the well-known problem describing the vibrational 
modes of a string with nonhomogeneous mass density. This fact lies at the core 
of our investigations of the inverse spectral problem. 



3.1 Liouville transformation to a finite interval 

A Liouville transformation is a change of both dependent and independent vari- 
ables in order to bring an ODE to its simplest form. The Liouville transforma- 
tion that we will use here was inspired by a similar transformation |2] relating 
the spectral problem for the Camassa-Holm equation to the string equation 
fiyyiv) = z 9(y) ( t>(y) with Dirichlet boundary conditions <^(±1) = 0. This self- 
adjoint spectral problem, which we briefly review in Appendix^ is fundamental 
in the inverse scattering approach to finding multi-peakon solutions of the CH 
equation As we will see, for the DP equation it is the non-selfadjoint 

"cubic string" that plays the corresponding role. 

To get an idea of how the transformation was found, consider the DP spectral 
problem (|2.31|l in the simplest possible case m(x) = 0, and leave aside the 
boundary conditions for a moment. Then the equation is just (d x — d^)if){x) — 
and the solution is 



ip(x) = Ae x + B + Ce~ x 



(e* + l) 2 



2e x 
2 



1 



y 



A 



2e x 



(i + yf 



B 



z B Ae x 
+ ~2 (e x + l) 2 

(l+y)(l-y) 



C 
2" 



c 



(1-1/) 



e x + 1 
2 



(3.1) 



where y = tanh(a;/2) 

y e (-1,1 



This change of variables maps x £ R to 



and the inverse transformation is x = log . The expression 
in brackets, call it <f){y), is a quadratic polynomial and thus satisfies —dy(f) = 0; 
the first order term in the operator d x — d x has been removed by the transfor- 
mation. 

In the general case, the following holds. 
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Theorem 3.1. Under the change of variables 

x , 2<j>{y) 

2' Ty ' 1-y 



y = tanh±, 4>{x) = f™, (3.2) 



i/ie DP spectral problem (|2.31|) is equivalent to the cubic string problem 

-<j>yyy(y) = z g(y) <f>(y) for ye (-1,1), 
0(-l)=^„(-l) = O, (3.3) 
0(1) = 0, 



g(y)=m(x). (3.4) 



/n i/ie discrete case, when m(x) — 2^2™ mi 6(x — Xi), equation i|3.4[l should be 
interpreted as 

n Sm- 
g(y) = ^2gi s (y - Vi), where y^tanhy, g t = % —^. (3.5) 

Proof. The first part is a straightforward computation with the chain rule, using 
|| = Y^y 1 ' ^ n particular, the equivalence of the boundary conditions follows 
from 

4>(±1)= Urn </>(*), 



b y (-l)= lim U a (a;) + i— ^V(a:)) 



For the discrete case, 8(x — Xi)dx = S(y — yi)dy, hence 

and the statement follows. □ 

By analogy with the ordinary discrete string, we will refer to the quantities 
gi as point masses at the positions yi, although physically speaking this termi- 
nology is perhaps not justified. We also define yo — — 1 and y n +i = 1, which is 
consistent with our convention that xq = — oo and x n +± = oo. 

Remark 3.2. For comparison, we remark that the Liouville transformation 
y = tanhx, ip(x) = (f>{y) /y/l — y 2 used in the CH case |2I maps %j){x) = Ae x + 
Be~ x S ker(9| - 1) to <j>(y) = A(l + y) + B(l - y) e ker(c^). 

Since the adjoint of l|3.3|) involves two boundary conditions at the right end- 
point and one at the left endpoint, l|3.3|l is not a sclfadjoint problem. Hence there 
is no a priori reason to suspect that the eigenvalues need to be real. However, 
this follows in the discrete case immediately from Theorems 13. II and 12.121 
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Theorem 3.3. The discrete cubic string problem l|3.3[) . 1|3.5[) . with all gi > 0, 
has n distinct positive eigenvalues 

< Ai < • • • < A„. 

In the general case, when the mass distribution g(y) is not discrete, there is 
an infinite sequence of eigenvalues 

< A : < A 2 < • • • , 

which we will prove in Section EOI using the beautiful theory of oscillatory kernels 
developed by Gantmacher and Krein. 

3.2 The discrete cubic string wavefunction <f)(y; z) 

Because of Theorem l3.ll all notions defined in the context of the peakon spectral 
problem (|2.31|1 have counterparts in the cubic string setting. The DP equation 
induces an isospectral deformation of the cubic string, but we mostly consider 
a fixed time t and suppress the time dependence in the notation. 

The DP wavefunction ip(x; z) of Definition 12.61 is mapped by the Liouville 
transformation (|3.2[1 to the cubic string wavefunction (j)(y; z), which by definition 
is the solution of the initial value problem 

<Pvvv(v) + z 9(y) ^{v) = for ye (-1,1), 

<H-i) = <M-i) = °> (3-6) 

Again, we are mainly interested in the discrete case 

n 

div) = ^9i&{y - Vi), 9i>0, 
1 

-1 = y < yi < ■ ■ ■ < y n < y n +i = +1- 

From H3.1fl it is clear that the peakon wavefunction, piecewise given by 

tp(x;z) = A k (z)e x + B k {z) + C k (z)e~ x , x € (x k ,x k+ i) 

according to (|2.20() . corresponds to the discrete cubic string wavefunction 

{l + y) 2 1 - y 2 (1-y) 2 
<j)(y;z) = A k (z) V B k {z) — — + C k {z) , y e (y k ,y k+1 ), 

(3.7) 

which is piecewise a quadratic polynomial in y. We will denote by 

h = Vk+i - Vk (3.8) 

the length of the fcth interval (y k , y k +i), for k = 0, . . . , n. The coefficients A kl 
B kl C k in H3.7JI are the same as before and are given by H2.26JI . In particular, 
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Aq = 1 and B = Co = 0, so that 4>(y; z) = i(l + y) 2 for — 1 < y < 1. However, 
we will not make much use of these coefficients from now on; instead we will 
keep track of 4>(y; z) through the values of the function and its derivatives at 
the points y k . We will use the vector notation 

/ <t>{y;z) 
®{y;z) = <l>y(y;z) 

sometimes omitting z to simplify the notation. 

By <|3.6fl , 4>yyy (y) is singular (proportional to the delta function) at the points 
yi, . . . , y n , and zero elsewhere. Thus we require <fi{y) and 4> y (y) to be continuous, 
and 4> yy {y) to be piecewise constant with jump discontinuities at y\, . . . , y n : 

<t>yy(yk+) - 4> yy {yk-) = -zg k (j)(y k ). (3.9) 

These jump conditions can be expressed as 




Hy k +) = G k (z)$(y k -), where G k {z) = | 1 | . (8.10) 
Further, 

<t>{y) = <t>(yk) + <t> v {yk) (y - yu) + ^yy(yk+) — ^ 

for y k < y < y k +\. Evaluating this expression and its derivatives at y = y k +\ 
tells us how $ propagates from one jump to the next: 

$(tfk+i-) = i**(!/fc+), where L k = 1 l k . (3.11) 

\0 1 J 

Clearly, the vectors Q(y k ±;z) will consist of polynomials in z (with coeffi- 
cients that depend on g^s and h's). The degrees are as follows, for k > 1: as 

z — > CO, 

$(y fc -;z) = L fc _ 1 G fe _ 1 (z)...L 1 G 1 (z)L = O^*" 1 ) (3.12) 

VV \0(z^)J 

and 

/0(z fe -!)\ 

$( ?A + ; z) = G k (z) $(y k -; z) = O^*" 1 ) . (3.13) 

V O(z^) / 

It is immediate that G k (z)~ 1 = G k (—z) and Lfc(^) -1 = L k (—l k ). Some slightly 
less obvious relations, which are crucial in the solution of the inverse spectral 
problem, concern the transposed inverses: 

(L^Y = JL k J and (G^z)- 1 )* = JG k {-z)J 1 (3.14) 
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where 




(3.15) 



3.3 Weyl functions 

The eigenvalues of the cubic string problem l|3.3|l are the zeros of the function 
4>(l;z), where <j>(y;z) is the cubic string wavefunction defined in the previous 
section as the solution of the initial value problem (|3.t)|) . The extended spectral 
data are encoded in a pair of Weyl functions. 

Definition 3.4. The Weyl functions for the cubic string are 

W(z) = *M Z{z) = (3.16) 
0(1; z) <P(l;z) 

where (f>(y; z) is the cubic string wavefunction. 

Theorem 3.5. In the discrete case g(y) = J^i 9k$(y — yk), 9k > 0, the Weyl 
functions admit the partial fraction decompositions 

m.i + ±-±--±* (3,7) 

Z Z — ' Z — Ah — ' Z — Afe 

fe=l fe=0 
Z(z) 1/2 c k c k ,„ 

fc=l K k=0 

where {Afe, bk, Ck}2— constitute the extended spectral data of the corresponding 
peakon spectral problem (see Theorem \3.1\ and Definition \2.1^\j . 

Proof. In the discrete case, 0(1; z) is an nth degree polynomial in z. The eigen- 
values Ai, . . . , A n , which are simple and positive by Theorem 13.31 are the zeros 
of <f)(l;z). This shows that W(z)/z and Z(z)jz have only simple poles, pre- 
cisely at z = = Ao and at each eigenvalue z — A&. To compute the residues 
at z = 0, we use that (/>(y;0) = \{1 + y) 2 , which gives b = = § and 

Co = ^({.q^ — \ m agreement with our definitions. 

In the rightmost interval y n < y < 1 we have, by (|3.7|) . 

<Kv,z) =A(z)^l + B(z) 1 -^ + Ciz)^ ~ »>' 



hence 



2 w 2 2 

W(z) _ <£„(!; z) _ 2A(z)-B(z) _ 1 



and 



z Z(f>(l;z) 2zA{z) z 2zA(z) 

Z(z) _ <j> yy (l; z) _ A{z)- B(z)+C(z) _ 1/2 C(z) - B(z) 

z ~ z4>(l;z) ~~ 2zA(z) ~ z 2z A(z) 

Comparison with (|2.33|) finishes the proof. □ 
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Corollary 3.6. The second Weyl function Z(z) is determined by the first Weyl 
function W(z) through equation l|2.36[) expressing Ck in terms of Ak 's and bk 's. 

This fact, which is absolutely fundamental to the solution of the inverse 
spectral problem, is far from obvious from the definition of the Weyl functions in 
terms of a given discrete cubic string. In fact, the proof of (|2.36|) in Section |2~7I 
was only made possible by knowing how the spectral data change under the 
peakon evolution. In other words, in order to prove Corollary 13 . 61 we needed to 
know about the isospectral deformation of the cubic string induced by the DP 
equation. 

3.4 Gantmacher— Krein theory and the cubic string 

Before we turn to the inverse spectral problem, we would like to show how 
the cubic string fits into the Gantmacher-Krein theory of oscillatory kernels 
|12|. We begin by collecting a few facts about three classes of matrices with 
nonnegative entries, see e.g QT] or [Ttj]. 

Definition 3.7. An nxn matrix A is called totally positive (totally nonnegative) 
if every minor of A is positive (nonnegative) . A totally nonnegative matrix A 
is called oscillatory if some power of it is totally positive. 

Oscillatory matrices can be thought of as being somewhere between totally 
nonnegative and totally positive matrices. The most pertinent property of these 
classes of matrices is the following. 

Theorem 3.8. All eigenvalues of a totally positive matrix are positive and of 
algebraic multiplicity one. All eigenvalues of a totally nonnegative matrix are 
nonnegative, but in general of arbitrary multiplicity. 

Corollary 3.9. All eigenvalues of an oscillatory matrix are positive and of 
algebraic multiplicity one. 

Definition 3.10. Suppose / is an open interval of R. A (necessarily positive) 
continuous function K(y,s), y, s G / is called oscillatory if for every choice of 
points yi < yi < • • • < y n € / the matrix [K(yi, yj)]ij = \ t ... t n is oscillatory. 

Consider the integral equation 



where the integral is taken in the sense of Stieltjes, with a a non-decreasing 
function on I. In the above, z plays the role of an eigenvalue and the sought 
solution <p is the corresponding eigcnfunction. The central result of the theory 
developed by Gantmacher and Krein is the following theorem. 

Theorem 3.11 (Gantmacher— Krein). If the kernel K of the integral equa- 
tion (|3.19|l is oscillatory, then: 




(3.19) 
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1. The eigenvalues are all positive and simple: < Ai < A2 < • • • . If the 
function o~(s) has only a finite number n of points of growth, then there 
are n eigenvalues. Otherwise, there are infinitely many eigenvalues. 

2. The eigenfunction </>i corresponding to the smallest eigenvalue \\ has no 
zeros in the interval I. 

3. For every j > 1, the eigenfunction cj>j corresponding to the jth eigenvalue 
Xj has exactly j — 1 nodal points in I (i.e., it has j — I zeros and changes 
its sign at each). 

The integral equation H3.19[l occurs most naturally in the context of bound- 
ary value problems for ordinary differential operators. 

Theorem 3.12 (Krein [14J). Let L{4>{y)) = ££ =0 l k (y) 0« (y), where l n (y) > 
0, and consider the boundary value problem 

L{4>) = 0, 

where p+ q = n > 2. Let G{y,s) be the Green function corresponding to this 
boundary value problem and let W(/i, fi,---, fk) denote the Wronskian of the 
k times differentiate functions fx, fi, ■ ■ ■ , fk- Then (— \) q G{y, s) is oscillatory 
if and only if there exist p solutions uii(y), ^2(2/), ■ ■ ■ , u p (y) of the boundary 
value problem 

L(<P) = 0, 

4(b) = cf>'(b) = ••■ = 0C«i-i)(6) = o 

and q solutions u> p +i(y), w p +2(y), ■ ■ ■ , w n (y) of the boundary value problem 

L(0) = 0, 

0(a) =0'(a) = ••• = <j) {p - 1) (a) = 

such that, for all y G I , 

u>i>0, W(wi,w 2 ) > 0, ... , W(L>i,UJ2,...,U n )>0. 

The cubic string equation (|3 . 3|) is a special case of this setup that corresponds 
to p = 2, q = 1, L = d 3 /dy 3 , and I = (—1, 1). The Green function of the cubic 
string is the solution of 

Gyyy(y,s) = 6 (j/ - s) , 

G(-l,s) = G v (-l,s) = 0, G(M)=0, 

namely 

(s- l) 2 (y + l) 2 

, -l<y<s<i, 

G{!, ^ = < (s + iny-lf (a 2 -l)(y 2 -l) ; -1 < s < y < 1. (3,20) 
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Lemma 3.13. Let G(y,s) be the Green function (|3.2U|) for the cubic string. 
Then K(y, s) = —G(y, s) is oscillatory. 

Proof. By direct computation one verifies that u>i(y) = 1 — y, ^(y) = — (1 — y) 2 , 
and u>z{y) = (1 + y) 2 satisfy the conditions in Theorem 13. 121 □ 

The cubic string problem Ij3.3|l is equivalent to the integral equation 

<Kv) = z J i K (y> s ) <^ s ) da ( s )' ( 3 - 21 ) 

where K(y, s) = —G(y, s) and da(s) — g(s) ds. In this formulation we can 
include any kind of (positive) mass distribution, by letting the non-decreasing 
function <x(s) equal the accumulated mass in the interval — 1 < y < s. Since K 
is oscillatory, Theorem 13.111 gives the following more general version of Theo- 
rem valid not only in the discrete case. 

Theorem 3.14. The spectrum of the cubic string problem (|3.21|l is positive and 
simple: < Ai < A2 < • • • . 

Theorem 13. 1 II also provides an alternative proof of the positivity of the bk's. 

Theorem 3.15. The residues bk in the Weyl function W(z) of the discrete 
cubic string are positive. (See equation B3.17JI . ) 



Proof. We have 



z=\ k zcf>(l;z) 2Utf h (l- Afc/Aj) 



where (f>(y, z) = 2 Jl™=i(l — Z /Afe) is the solution to the cubic string initial value 
problem with </>(— l;z) = 4> y {— I; z) = and <p yy (l;z) = 1. In other words, 
(f>(y; Afe) is the fcth eigenfunction of the cubic string boundary value problem, 
and as such it has k — 1 nodal points according to Theorem l3.11l It follows that 
(— l) k (j)y(l; Afc) > 0. Suppose equality holds. Then cj)(y; Afe) is simultaneously 
an eigenfunction of our usual cubic string and of the mirrored cubic string 
(j> = 1, q = 2) with boundary conditions 0(1; z) — (j> y (l] z) = z) = 0. This 

is a contradiction, since it follows immediately from Theorems l3.11l and l3.12l that 
the mirrored cubic string has negative eigenvalues. Hence, (— l) k (f>y(l; Afc) > 0. 
And since precisely k — 1 factors in the denominator of 13.22|) are negative, it 
follows that bk > 0. □ 



4 Inverse problem for the discrete cubic string 

Our aim in this section is to give an explicit solution of the inverse spectral prob- 
lem for the discrete cubic string (|3.3|l . (|3.5|l . This will also provide us with what 
is needed to complete the proof of Theorem 12.221 regarding the corresponding 
inverse spectral problem for DP peakons. 
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The inverse problem for the ordinary (not cubic) discrete string plays a 
fundamental role in solving the peakon problem for the Camassa-Holm equation 
For the reader's convenience, we sketch the Stieltjes-Krein solution to 
that inverse problem in Appendix [S] with emphasis on the structures common 
to both inverse problems. 



z 

i=0 



4.1 General setup 

We begin by stating the problem precisely. 

Definition 4.1 (Inverse problem). Consider a discrete cubic string specified 
by a positive measure g(y) = 9* ^Viiv) wnose support satisfies — 1 < yi < 

?/2 < ■ • ■ < Vn < 1 ■ The inverse spectral problem is to determine the measure 
g{y) given the Weyl function W(z) of the cubic string: 

W{z) _ bi _ fdn{\) 

where 

n 
i=0 

fro = 1, bk > 0, and = Ao < Ai < • • • < A„, 

Recall that by Definition 13.41 the Weyl functions are 

<t>{l;z) 4>{l;z) 
where, according to ^FWfy , $FTfy . 

*(l;z)= <l> y (l\z) =L n G n (z)L n _ 1 G n - 1 (z)---L 1 G 1 {z)L o \0 
WM)/ V 1 / 

Also remember that Z is determined by W ( Corollary 13. 6|) . 

The analogy with the case of Stieltjes continued fraction (see Appendix lA"|l 
suggests that the 3x3 matrices 

a^(z)=L n 
a<V(z)=L n G n (z) 
a<V(z) =L n G n (z)L n _ 1 
a,W(z)=L n G n (z) L n ^G n ^{z) 
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might contain convergent-like quantities from which one can construct succes- 
sively better rational approximations to W(z) and Z(z). Likewise, from the 
vectors 



$(yo+;z)= (0,0,1)' 
$(lft-;z)= L (0,0,1)* 

$(yi+;z)= Gx{z)L (0,0,1)* 

$(y 2 -;z) = LiGi(z)L (0,0,1)* 

$(y 2 +; z) = G 2 (z) L x d(z) i (0, 0, 1)* 



one should be able to form remainder-like quantities. This is indeed the case, 
as we shall soon show. 

In contrast to the classical Stieltjes theory, it is not necessary to consider 
both the even and the odd case. This is because the matrices a^ 2k \z) and 
a( 2fe_1 )(z) differ only in the first column, and it turns out that the information 
that we want to extract can be found in their common second column. We 
choose here to study the even case, where the starting point is the relation 

$(l;z) = a( 2fc )(z)$(jM_ fc+1 -;z). (4.1) 

We will for the most part regard k 6 {1, . . . , n} as fixed (but arbitrary), and 
denote the matrix a^ 2k ^ (z) simply by a(z) or a. We collect some basic facts about 
it in the following lemma. As in Section lTHl it is convenient to let k' = n + 1 — k. 

Lemma 4.2. The matrix 

a{z)=a^ k \z)=L n G n {z)L n ^G n ^{z)---L k ,G kl {z) (4.2) 

satisfies deta(z) = 1. The entries in its first column are polynomials in z of 
degree k, while the remaining entries are polynomials of degree k — 1. In the 
second column we note in particular that 

M*) = h- ( ft ^f] zk ~ X + ■ ■ ■ + ( E l ) ( 4 - 3 ) 

\i=fe'+l / \i=k' J 

and 

022(0) = 1, a 32 (0) = 0. (4.4) 
All 2x2 minors of a(z) are polynomials of degree at most k. 

Proof. The matrices Gi(z) and Li have determinant one, hence so does a(z). 
Since Gi{z) is linear in z, it is clear that a(z) is a matrix polynomial in z 
of degree k. However, the second and third columns in a(z) are unaffected 
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by the last factor Gk'(z), and are consequently only of degree k — 1. To ex- 
tract the leading coefficient, write Gi{z) = I + z (— <?i)(0, 0, 1)*(1, 0, 0) and use 
(1,0,0)^(0,0,1)' = if/2. The constant term is 




o(0) = L n £„-i • • • = | o ! £^1, 

1 

Finally, because a(z) has determinant one, its adjoint equals a(z) , which is 
of degree k since Gi(z)^ 1 — Gi(—z), thus implying that all 2 x 2 minors of a(z) 
are polynomials of degree at most k. □ 

Returning now to equation (|4.1|l . we divide it by 0(1; z) and let $(?/£/ — ; z) = 
(Pk,qk,r k y to obtain 

/p*(*)\ 

a< 2fe )(z) . (4.5) 

Suppressing the dependence on z and k in order to simplify the notation, we 
find 

W a 2 ip + a 22 q + a 23 r a 2 i R + a 22 + a 23 R , . 

— = = 7s, (4.6) 

1 a n p+ a 12 q + a 13 r a\i R + a 12 + a 13 R 

where R = p/q and R = r/q. Similarly, 

Z_ _ a 3 i R + a 3 2 + a 3 3 R 
1 an R + ai2 + ai3 # 

The quantities i? and R, or more precisely 

^ = ^ = -^4 and % h = r ±= ^'->*\ (4.8) 

are the analogs of remainders referred to earlier. By formally setting the re- 
mainders equal to zero in (£333) and g2J 

we get the analogs of convergents, 
W w a 22 /a\ 2 and Z ss 032/012, with the following order of approximation: 



(4.7) 



Theorem 4.3. As z — > 00, 



a 12 ( Z ) 

a 12 (zj 
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Proof. It follows from (|4.6|l that 

012021 — 011022 ai2»23 — 0130-22 S 

ii H rt 

™. _ 022 _ Ql2 012 

Ol2 On i? + Oi2 + Oi 3 i? 

By l|3.12|) . p, q, and r are all of degree k' — 1, hence i? = p/q = 0(1) and 
i? = r/q = 0(1). From this, and from Lemma f4. 21 we see that the right-hand 
side equals 

°^o(i)+°£Lo(i) 



0{z k )0{\) + 0(z fc -!) + o^-^cki) 
The proof for Z is similar. □ 
Theorem 4.4. Let 

W*(z) = -W(-z), Z*{z) = Z(-z). (4.10) 

Then, with a = a 

Z* ai 2 + W*a 2 2+a 32 = \J^J , asz-xx. (4.11) 
Proof. From 14.5( 1 we have 

a(z)- 1 I 1 : : —— <!'(////■ ■::) : < M TT 



Z(z) 
impli 

(Note that J = J -1 = J'.) Hence 



0(1; z) 

Now recall (|3.14H . which implies that (a(z) -1 )* = Ja(—z)J, where J = [o -io 



Ja(-z) 4 J (V(z)| =0(-V) 



which upon changing z to -z and multiplying both sides by J yields 

Z(-z) 

a- 
i 

the second row of which is precisely 14.11|l . □ 



a{zf I -W(-z) \=0[\), 
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4.2 An approximation problem 



Lemma |4~21 Theorem 14. 31 and Theorem 14 .41 imply that the entries in the second 
column of a = a^ 2k ^(z), 

(Q,P,P) = (a 12l a 2 2,a 32 ), 
satisfy the following approximation problem: 

Definition 4.5 (Pade-like approximation problem). Let the functions W 
and Z be given by 

Z L> Z-X k ^ Z3 + 1 ' 

(4.12) 



z ( z ) _ \ - ^fc = \ - ij 

z ^ z-X k ^ zi+ l ' 

fc=0 j=0 



where 



& = E^ A i 7,- = E c ^ ( 4 - 13 ) 

fc=0 fc=0 

and the c^'s depend on bkS and A^'s as in Corollary 12. 161 For a given integer 
1 < k < n, we seek three polynomials (Q, P, P) of degree k — 1 satisfying the 
following conditions: 

1. (Approximation) 



2. (Symmetry) 



Z*<3 + M / *P + P = 0( ^t:) (z^oo), 



where W^(z) = -W(-z) and Z*(z) = Z{-z). 
3. (Normalization at z = 0) 

P(0) = 1, P(0) = 0. 



Remark 4.6. The approximation condition is similar to classical Pade-Hermite 
approximation |17|. where one approximates two (independently chosen) func- 
tions to order 0{1/ z 2 ^ 1 ^) using a common denominator of degree k — 1. Since 
in our case we only have 0(l/z' c ~ 1 ), there is clearly not enough information in 
the approximation condition alone to uniquely determine the denominator Q{z). 
However, our W and Z are not independent of each other, and the additional 
symmetry condition will eventually lead to a unique Q(z). 
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If we introduce spectral measures 



M = E &fe( ^ fc and v = ^c k Sx k , (4-14) 

k=0 k=0 



then 



(3j = 6feAfe = / x j dn(x) and = CfeA{, = / x j dv(x) (4.15) 

fc=0 fc=0 

are the moments of /i and is, respectively. By Corollary 12. 161 v is determined 
by fi. In fact, l|2.36l) can be expressed as 

f(x) dv{x) = [[ X -^± d^x) dn(y) = [[ d»{x) dfi(y). (4.16) 

In particular, 

lo = II : —~ dn(x) dfi{y) = 1 1 f— - dfi{x) dn{y). (4.17) 
JJ x + y JJ x + y 

Lemma 4.7. The relation 14.16f) between /i and v implies that W and Z satisfy 
the constraint 

Z(z) + Z*(z) + W*{z)W{z) = 0. 
Proof. Expressing W and Z in terms of the spectral measures we obtain 

z 

z + y 



x + y J J x + y 



x J+1 f f y j+1 



W(z) = I ^— d»(x), W*(z) = - / — dfx(y), 



and 



f z f f zx 1 

Z(z) = I du(x) = // ■ — dfi(x)dfi(y), 

J z-x JJ z~xx+y 

/z if zy 1 
dv(x) =11 du(x) du(y), 
z+x J J z+y x+y 

from which the identity in question follows immediately. □ 

For the purposes of this approximation problem, we can consider all functions 
as being formal Laurent series with finitely many positive powers. Let II± and 
II^ denote the following projection operators on finite dimensional subspaces 
(we regard k as fixed, and omit dependence on k in the notation): 

, N v fe-i , N s fc-l 

n +( E = E a ^ n +( E = E a * z4 > 

( N \ 
n_f a i z ) = E a%zl - 

X »=-(k-l) 



i — — oo 
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Let Mf denote multiplication by f(z), and define the truncated Hankcl and 
Toeplitz operators 



H f = n_ o Mf o n° : span{z°, . . . , z^ 1 } -> span{z _fe+1 , z }, 

z k ' x } -> spanjz 1 , . . . , z^ 1 }. 



(4.18) 



Tf =11+0 Mf o 11+ : span{z°, . 
Moreover, let 

i : spanjz 1 , . . . , z k ~ 1 } — > span{z°, z 1 , . . . , z fc_1 } 
be the natural inclusion. 

Theorem 4.8. TTie polynomials Q(z), P(z), P(z) solve the approximation prob- 
lem of Definition ^. 5\ if and only if 



P = 1 + U+WQ, P = Il + ZQ, 
[H z * + H w * o i o T w ] Q = -IL.W* 

Equation (|4.2U|) is equivalent to the linear system 



(4.19) 
(4.20) 



( <5oo 
<5io 



<5oi 
Su 



$o,k-i \ 



5k-l,h-\J \1k-l) 



(4.21) 



/or i/ie unknown coefficients in Q(z) — ^2 k = Q ft- 2 *' where 



Sab 



x a+1 y b 
x + y 



dfj,(x)dfj,(y). 



(4.22) 



Proof. According to the approximation condition, the functions WQ — P and 
ZQ — P contain no positive powers of z, hence their 11+ projections are zero. 
This, together with the normalization conditions, is equivalent to l|4.19[) (which 
shows that P and P are uniquely determined by Q). Inserting (|4.19(l into the 
symmetry condition yields 



z* q + w* (i + n+wQ) + n+zQ = o (jA 



(4.23) 



The coefficients of the positive powers z 1 , 



.,/,■-! 



appearing in (|4.23() vanish 



identically regardless of Q; indeed, using (11+ — I)Q = ~qo, Tl+W* — 0, and 
Lemma 14.71 we obtain 



n+ [z* q + w* (i + n+ wq) + n+zQj = n+ [z + z* + w*n+ wj q 
= ti + (z + z* + w*(n + -I + T)W}Q = n+(z + Z* + W*W^Q = 0. 
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On the other hand, the vanishing of the coefficients of z ^ l \ . . . , z° in H4.23[) 
is equivalent to the following condition on Q: 

U (z* Q + W* (1 + U+WQ) + U+ZQ 



0. 



(4.24) 



which is just (|4.20() . 

To show equivalence of (|4.2U|I and (|4.21|) we write out l|4.2U|l in terms of 
matrices with respect to the standard bases ordered as in (|4.18|l . Then Q is 
represented by the column vector (go, ■ ■ ■ , <Zn— l) j while 



(-l) a+b la 



+b 



0<a<k-l 
0<6<fc-l 



/ 7o -7i 72 
-7i 72 
72 



V 



(-1) 



a ■ Hi + 1 



Pa, 



0<a<fc-l 
Km<fc-1 



(Pi -02 03 
-02 03 
03 

V : 



Tw = 







b—m 



l<m<k-l 
0<6<fe-l 



(0 00 01 02 

o o A A 

p 



-I1_W* = (P , -01,02, ■ • ■ , (-l f~ L k -l) - 

If we multiply out the matrices and change the sign of every second row, the 
system takes the form 

fc-l b 

[(-l) b 7a+h - (-l) m 0u+m0b- m ] q b = 0a, < a < k - 1. (4.25) 

fc=0 m=l 

(When 6 = 0, the inner sum is empty.) Upon using the spectral representations 
H4.15fl and (|4.17|l of P's and 7's we arrive at the claimed integral representation 
H4.22fl of the (a, b) entry 8 a b of the matrix of coefficients in l|4.25|l : 



Sab = 



(-1) 6 



r a+b+l 



x + y 



- J2(-l) m x a+m y b - ,n 



dn(x) dfj,(y) 



x a+l y b 

x + y 



dfj,(x) dn{y). 



□ 



To complete the picture, it remains to show that the system (|4.21(l has 
nonzero determinant, so that the solution to the approximation problem exists 
and is unique. This will be dealt with in the following sections. 
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Remark 4.9. The choice of convergents and remainders is not unique. For 
example, one could take the entries in the first column of a^ 2 ^ (z) as convergents, 
and instead of (|4.6[) write 

q r 
&21 + «22 - + dl3 - 



a ll + a 12 h Ql3 — 

P P 

where now q/p and r/p play the role of remainders, and the order of approx- 
imation can be shown to be 0(1/ z k ); similarly with Z and equation Q4.7|) . In 
this way one obtains an approximation problem which is similar to, but slightly 
different from, the one in Definition 14.51 When looking at the first column, the 
odd matrices a( 2fe_1 ) are different from the even ones, and this case provides 
yet another approximation problem of similar type. However, we are convinced 
that the problem we have chosen to focus on here is the simplest route to the 
solution of the inverse problem. 




4.3 Determinants 

From the expression 14.22JI for 5 a b as a double integral, it is clear that the 
determinants that appear when one tries to solve the linear system (|4.21|l with 
Cramer's rule are somewhat reminiscent of the classical fcx fc Hankel determinant 
of moments, 



0o 
0i 



01 

02 



0k-i 
0k 



0k-l 0k 



0: 



■2k-2 



1 
fc! 



A{x) 2 dfi k (x), 



(4.26) 



R 1 



where A(x) = A(jci, . 



, Xk) 



n 



(For a proof of Heine's formula 



H4.26fl . see for instance 10., page 27] or [SJ Proposition 3.8].) 

The following lemma contains similar integral formulas for the present deter- 
minants. These formulas are valid for any measure, although our main interest 
here is of course in the discrete case when \i = X^T=o 



Lemma 4.10. Suppose fi is a measure on R such that the integrals 



f3 a = J x a dfi(x) and S a b 
are finite. Let 



x a+l y b 



dfi(x) dn(y) 



(4.27) 



Uk 



1 

fc! 

Vk = h 



R 



A(.t) 2 

A(x) 2 
it* T ( x ) 



dn k {x), 



X\X2 ■ ■ ■ Xk d/i k (x), 



(4.28) 
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wher 



A(x) = A(xi, . . . ,x k ) = W[xi - Xj), 

r(x) = r(xi, ...,x k ) = Y\_(xi + Xj). 



(4.29) 



i<3 



(When k = or 1, we let A(x) = T(x) = 1. Also, uq = vq = I.) Then for k > 1 
the following k x k determinant formulas hold: 



and 





5oo 


5ox 


$02 


5 ,fc- 






Su 


5x2 


0~l,k- 


D k := 


^20 


$21 


522 


5 2 ,fc- 




5fe-X,0 


dk-i,i 


Sk-1,2 ■ 


5fe-x,fe 




A) 


5oo 




5(),fc-2 




/3x 


5xo 


Sn .. 


5x,fc-2 




ft 


^20 


$21 


52,fc-2 




Pk-l 


Sk-i,o 


Sk-1,1 ■ ■ 


5fe-X,fe-2 






Soi 


5()2 


5o,fc-x 




Pi 


Su 


5x2 


5x,fc-x 




fa 


$21 


522 


52,fc-X 




Pk-l 


Sk-i,i 


5fc-X,2 • ■ 


5fc_i i fc_ 



2 fc 



(4.30) 



o/c — X : 



(4.31) 



UkVk-l 

2 k~i 



(4.32) 



Proof. We show only l|4.31|l ; the formulas (|4.3U|) and (|4.32|) are proved similarly. 
Notationally it is slightly easier to work with the (fc + 1) x (fc + 1) determinant 
D £ , so what we will actually prove is that 



D'k+i 



Uk u k +i 
2 k 



(4.33) 



Note first that our definition of A(a;) is such that the Vandermonde deter- 
minant equals 

det(^" 1 )i J= i ! ... ! fe = Y\_( x i - x j) = A(x fc , ...,Xi) = (-l)( 2 )A(a;i, . . .,x k ). 

In the final result the sign convention does not matter since only A(x) 2 appears 
there, but we need to keep track of the correct signs during the proof. As regards 
r(x), the order of the variables is of course immaterial. 

Like in the proof of Heine's formula, we use different dummy variables in 
each column, in order to be able to pull integrals outside of the determinant 
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by multilincarity. We will need 2k + 1 variables (xo,Xi, . . . ,X2k), &s follows. 
In the integrals (3 a in the first column we call the integration variable xo, in 
the double integrals <5 a o in the second column we use x\ and x 2 , in the double 
integrals b a \ in the third column we use x 3 and X4, etc., with even-numbered 
variables replacing x and odd-numbered variables replacing y. Then we can pull 
out all integral signs, and also the factor X26-1 x 2 b / (2:26-1 + X2b) from column 
6+1, for b = l,...,k. What remains is a Vandermonde determinant in the 
even-numbered variables, so 

„, f {x\x\xl- ■ ■ X k 2k ^ l ){x 2 X A - ■ ■ X 2k ) A(x 2 fc, . . . ,X 4 ,X 2 ,X ) 

D k+i = / 7 w ; 7 n cfyr* +1 (x). 

J R 2k+i (Xi + X 2 )(x 3 + Xi) ■ ■ ■ {X 2 k-1 + X 2 k) 

We now introduce an auxiliary variable X-\ = in order to write the last 
expression in a more symmetric form: 

D , I i x l x l x l ■ ■ ■ X2kl 1 ){x Q X 2 X i ■ ■ ■ X2k) A(x 2 fc, ■ ■ ■ , £4, X 2 , Xp) ^ 2fc+l fa) 

k+1 7 R 2fc+i (x_i +X )(X1 + X 2 )(x 3 + Xi) ■ ■ ■ {X2k-1 +X2k) 

The value of the original determinant is independent of the order in which 
we number the dummy variables, so the above expression is invariant under per- 
mutations of the variables. In particular, symmetrizing over the even-numbered 
variables doesn't change the value. That is, if Seven denotes the set of permu- 
tations of {0, 2,4,..., 2k}, then 

D' k+l = . 1 - / (xlxl---X 2 C kl 1 ){x X2X 4 ---X2k)A(X2k,---,X4,X2,X ) 

(fc + 1)! J R 2fc+l 

J y ( s § n7r ) \ d^ k+1 { x ) 

L c q (x-i+x 7r{0) )(x 1 +x 7r{2) )(x 3 +x 7r(i) )---(x 2 k-i+x n{2k) ) 

\7TtOeven / 

The sum in parentheses is a Cauchy determinant, which by the well-known 
formula equals 

AQ_i, xi, x 3 , x 2 k-i)A(x ,x 2 ,x 4 , ■ ■ ■ , x 2k ) 

nLo n}=o foi-i +X2j) 

_ (-xi)(-x 3 ) ■ ■ ■ (-x 2fc -i)A(xi,x 3 , . . . ,x 2 fc-i)A(xo,x 2 ,X4,- ■ ■ ,x 2 k) 
(x x 2 x 4 • ■■X2k) rij=i Ilj=o( x 2i-i + x 2j ) 

where the second line is a result of splitting off terms involving x_i. It follows 
that 

D *+1 = 7FTW I {-l) {ktl) + \x l X Z ---X2 k - 1 ){xlx\---X k 2 -\) 

(It + 1J! J R 2fe+l 

x A(xi,x 3 , . . . ,x 2fc _i)A(x ,x 2 , . . . ,x 2fc ) 2 dfj 2k+i^ 

ULlU k 3 =o(x2t-l + X 2 j) 
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Symmetrizing this over the odd-numbered variables we get yet another Vander- 
monde determinant (— 1)( 2 ) A(xi, X3, . . . , x%k-i)i hence 



and then symmetrizing over all variables gives 



A(xi, x 3 , . . . , x 2 k-i) 2 A(x , x 2 , ■ ■ ■ , x 2 k) 2 d(J 2k+i/ x \ 



D ' k+1 ~ k!(k + l)\(2k+l)\J R2k J S (^(D^O) •••^(2fc-i)) 



A(x 7r ( 1 ), Xn@y, Xn^k-l)) 2 A(x ff (p), a; ff ( 2 ), ■ ■ ■ , X 7r (2k)) 2 \ , 2k+l( \ 

1 11^=0(^(2.-1) + 2^(23)) / 

If (T G Safc+i leaves the sets {0,2,..., 2k} and {1,3, ... ,2k — 1} invariant, then 
the terms corresponding to tt and to tt o cr are equal, so each term actually 
appears k\(k + 1)1 times in the sum above. Using the notation (}°C^) for the set 
of fc-element subsets {«!<•••< ik} of {0, 1, . . . ,n}, as well as other notation 
introduced earlier in Definition 12. 181 we can write D' k+1 in a less redundant way 
by picking out a "sorted representative" from each such group of equal terms: 



(2k + iy. j R 2 k+1 \f^ iyj 

/ ^(^AjA^W+V), 
j R 2fc+i 1 iuj J 



(2fe + l)! y R2t+ i r 7 



where the sum runs over all ( 2k ^ 1 ) ways of partitioning {0,1,..., 2k} = I U J 

into disjoint sets J € an d ^ G ('fe+i')- Now we claim that the symmetric 

polynomial given by the sum in parentheses satisfies the identity 

with / and J running over the same sets as before. The general strategy for 
proving identities of this type is outlined in Appendix [Bj which also contains 
other identities needed in the proof of the remaining two statements (|4.30ll and 



47 



14.32fl of the present theorem. Granted the claim, we obtain 

/ R 2fc + 1 



D' k+1 = 



(2fc + l)! 

2~ k k\(k- 
(2fc+l)! 



E 

7, J V 



V fe+1 (z) 



(fc + 1)! 



A 2 

R<=+1 1 J 



-J 



f2fc+ 
I k 



K^2 U k Uk+1 
I I, J 



Uk u k +i 

2 k ■ 



which proves H4.33JI . 



□ 



Here we need yet some more notation. Let V—\ = 0, Uq = Vq = 1, and for 
k > 1: 



(A/)- 



In, 



E 

/6 (:o,»]j 



Tr 



Xibi, 



(4.34) 



In other words, Uk and Vfc are symmetric functions of the same form as Uk 
and 14. of Definition 12. 181 but involving the additional variables Ao and b . For 
future use we also define, for k > 0, 



U k 

Uk+l 



Vk-i 



(4.35) 



Lemma 4.11. When [i — X)™=o^ i( -* A ;> ^he integrals Uk and Vk in 14.28|l reduce 
to the sums Uk and Vk, respectively. Hence, in this case, 

> 2 



(UkY 



D' k = 



UkU k - 



D' k ' = 



UkV k - 



2k ' * 2 k ~ 1 ' 2 fe ~ 1 ' 

Proof. This follows from the definitions and Lemma f4. 101 

Recall that in the context of the discrete cubic string we have Ao 
6o = 1 by definition. 

Lemma 4.12. When Ao = and bo = 1, 

U k = U k + 14-i, V k = V k , W k = W k . 



(4.36) 
□ 

and 
(4.37) 



Proof. In the sum defining Uk, the terms for which 0^1 add up to Uk, while 
if / = {0} U J with J e (fi^), then (A/) 2 /Tj = A, 7 (Aj) 2 /T, 7 and 6 7 = bj, so 

those terms add up to Vk-i- For Vk, the terms with £ / vanish, leaving only 
terms which add up to Vk ■ Finally, 



W k = 



U k + Vk-i 
Uk+i + v k 



Vk-i 

v k 



u k 

Uk+i 



Vk-i 

v k 



= w k . 



□ 
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4.4 Solution of the inverse problem 

Now we merely have to put the pieces together. 

Theorem 4.13. The approximation problem of Definition has a unique 
solution: 

1. The coefficients of the polynomial Q(z) = ^2i=o qiZ 1 are uniquely deter- 
mined by the nonsingular linear system H4.21|l . In particular, the highest 
and lowest coefficients of Q(z) are given by the following ratios of deter- 
minants: 

(-l) fc "V-i = qo = (4.38) 

2. The polynomials P(z) and P{z) are uniquely determined by Q(z) as 

P=1 + T W Q, P = T Z Q. (4.39) 

Proof. The determinant of is D k = (U k ) 2 /2 k = (U k + V k -i) 2 /2 k , by 

the lemmas in the previous section. This is clearly positive since U k > and 
Vfc-i > 0, so the system is nonsingular. Solving 14.21(1 for qo and q k -\ using 
Cramer's rule, one obtains the k x k determinants in Lemma [4.101 the factor 
(— l) fe_1 appearing when moving the /Vs from the last column to the first in 
the case of q k -i- Equation H4.39|l is just a reformulation of (|4.19|) . □ 

Expressing the result in terms of C/'s and V^'s using H4.36|l yields 
Corollary 4.14. 

(-1) Qk-i = * = -~— • (4.40) 

Uk <J k 

Remark 4.15. Note that there is a factor of U k which cancels in the quotients 
H4.38fl . As should be clear from the proof of Lemma T4. 101 this is a considerable 
complication compared to the classical case (Stieltjes continued fractions, Pade 
approximation, ordinary discrete string, Camassa-Holm peakons). 

Returning to the discrete cubic string, recall that (Q,P,P) = (012,^221^32) 
is a solution of our approximation problem. Comparing equation (|4.3|) for 012(2) 
to equation l|4.40|l for Q(z), we see that 




2V k . 



Uk + Vk-i 

2{Uk-i + V k - 2 ) 
Uk + Vk-i ' 



(4.41) 
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Given the quantities in 14.41f) one can easily compute the original string variables 
from 

1 V* I ; 2 Ofc+l 

Vk' = 1 - > J k, h' = V(k-iy ~ Vk>, gy = i r~FT> 

r£ i(k+i)'lk>Mk 

with the following result: 

Theorem 4.16. The solution to the inverse spectral problem for the discrete 
cubic string in terms of symmetric functions of A 's and b 's is 

1 2V k -i Uk-Vk-i 1 . 

Vk ' = = Uk~+Vk~i ( fc = 1 —^ 

2W k -i 2W k -! 
L> = — — ~ -, w r (k = 1, . . . , n + 1), (4.42) 

Uk-iUk {Uk-i + Vk-2)(Uk + V k -i) K ' ' h 1 ' 

(Uk) 4 _ (Uk + Vk-i) 4 , 
9k '~2Wk- 1 W k ~ 2Wk-iW k (k-l,...,n), 

with Uk, Vk, and W k as in Definition \2.1t<\ 

Remark 4.17. The formula for yk> is valid also for k = and k = n + 1; since 
Vq = U n+ i = 0, it yields y n +i = +1 and y = —1, respectively, as it should. 

Now we can finally prove the remaining part of Theorem 12.221 namely the 
explicit formulas for the inverse spectral problem in terms of peakon variables 
Xk and mfc. 

Lemma 4.18. 

x k <=\og- — , m k > = (k = l,...,n). (4.43) 

Vfc-l WkWk-l 

Proof. By l|3.5|l the string variables can be expressed in terms of peakon variables 
as 

x k , = logl±^, my = \g k ,{l - yl)\ {AAA) 

which immediately yields (|4.43l) upon inserting 14.4211 . □ 

So far the solution of the inverse problem has been obtained under the as- 
sumption that the spectral measure is coming from a cubic string. Theorem l2.22l 
implies however that there is a bijection between discrete cubic strings and spec- 
tral measures \x. 

Theorem 4.19. For every discrete cubic string specified by a positive measure 
9(y) = YJi=i9iby,(y) whose support satisfies -1 < j/i < y 2 < • ■ • < Un < 1 
there exists a unique positive measure 



^biS Xi , (4.45) 



i=0 
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bo = 1, bk > 0, = Aq < Ai < • • • < \ n , such that the Weyl function W(z) of 
the cubic string satisfies 



Conversely, given a measure [i as in (I4.45|) there exists a unique cubic string 
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A Appendix: The discrete string 

Here we give a brief account of the solution of the inverse spectral problem 
for a discrete string with Dirichlct boundary condition, for comparison to the 
discrete cubic string treated in this paper. For more details, we refer the reader 
to [2], and to Stieltjes' famous Recherches sur les fraction continues, especially 
the introduction and Chapter 3 Vol. II]. It was M. G. Krein who pointed 
out the mechanical interpretation of Stieltjes' work and generalized it in his 
spectral theory of the general (not necessarily discrete) inhomogenous string 
jl.'jj : see also Supplement II in or the comprehensive treatment in 0. 

Small vibrations u(y,t) of a string with mass density g(y) are described 
by the wave equation u yy = g{y)u u - If the ends of the string are attached 
at y = ±1, separation of variables u(y,t) — 0(y)r(i) results in the following 
equation for 4>{y): 



This is a selfadjoint spectral problem with simple, negative (if g > 0) eigenvalues 
z = \k = where tOk is the frequency of the fcth harmonic. The discrete 

string consisting of point masses gi at positions —1 < y\ < ■ ■ ■ < y n < +1 
corresponds to g(y) = giS(y—y,i) being a discrete measure. In this case, there 
are exactly n eigenvalues, and the eigenfunctions 4>{y) are piecewise linear (as 
opposed to the well known sinusoidal shape in the homogeneous case g(y) = 1). 
The piecewise constant slope 4> y {y) jumps by z gi 4>(yi) at each point yi. 




defined by Theorem ^. lb\ 



^yyiy) 



) = z g(y) <t>{y) for ye (-1,1) 

<M-i) = o, 0(1) = o. 



(A.l) 
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Figure 1: Discrete Dirichlet string with three masses 



The spectral problem for the discrete string can be conveniently reformulated 
using the transition matrix associated with the initial value problem posed at 
the left end y = —1. One starts at <p(— 1) = with initial slope <f> y {— 1) = 1, 
and propagates the solution to the right end point, which results in 

■KM) \. _ T r> (~\ t n i„\ j _ f° 
cj>y{l;z) 



L n G n (z) ■■■L 1 G 1 (z)L G\ , (A.2) 



where 

i)' <*<*>- ut 

(Like for the cubic string, we set Ik = Vk+i — Uk, and also k' = n + 1 — fe.) The 
eigenvalues are the zeros of the nth degree polynomial 0(1; z), and consequently 
they are the poles of the Weyl function of the problem, defined by 

W{z) = Miif). (A.3) 



With 



Qo Qi \ _. (i) 



L n — ( ) — ■ a 



L.,G n (z)=[ Q pl (A.4) 



L n G n (z)L n _ 1= ^ t) =: ° (3) ' 



and so on, equation (|A.2f) is just a matrix version of the standard recurrence for 
the convergents 



_ _Fb 1_ _ J\ 1 P 2 1 P; 

1 Qo In Ql , 1 



zg 



2 1 ^3 
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of a continued fraction expansion of Stieltjes type |18| . 
ln + 



W(z) = = £±1. (A.5) 

1 Wn+l 



zg n 



1 



ln-1 + 



Z92 



1 

h + 



1 



The remainders 



1 1 1 

Rq = — , i?i — — , i?2 — " , 

'0 1.1 

zgi + j- h H r 

'o 1 
251 + -r 

in the continued fraction are given by 

R 2k _ 1 = !!i, R 2k = ^, where ( q A = ( 5 *> ). (A.6) 

Pk q k +i \PkJ \4>y{yk-\z)J 1 ' 

Note that the convergents (Pj , Qj ) are obtained by formally setting the remain- 
der Rj to zero in Write 

fc=0 i=0 

where Ao = 0, ao = 1/2, and where 

n « 
fc=0 

is the jth moment of the spectral measure 

n 

f i = Y,akS-x k . (A.9) 



fe=0 



The convergents are Pade approximants of W: 

u// , -P 2 fc-i(z) n ( 1 



hk-t(z) \z 2 
degP 2 fe_i = fc-l, degQ 2 fc-i = ft - 1, P 2 fc-i(0) = l, 
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in the odd case, and 



Pik(z) . / 1 



w{z) - j^y- = ohi 

degP 2 fe = fc-l, degQ 2 k = k, Q 2 fe(0) = 1, 

in the even case. These conditions uniquely determine the coefficients of Qj(z) 
in terms of certain Hankel determinants AjL involving the moments Aj appearing 
in the Laurent series of W{z). Since 

Q2k-i(z) = (g n l n 9n-i ln-i ■ ■ ■gk')z k ~ 1 + ■ • • , 

Qlk{ z ) — (<?n In <?n-l ' ' ' fffc' ^fc') 2 ^ + . . . , 

we can express the string data as ratios of coefficients of Qj's, which themselves 
are ratios of Hankel determinants, resulting in altogether four determinants. In 
the notation of [2], 

, _ ( A fc-i) 2 B (A° fc ) 2 

These Hankel determinants can be evaluated in terms of A^'s and a^s using 
Heine's formula l|4.2(j[l . thereby recovering the string data {gi,U} in terms of 
symmetric functions of the spectral data {A^., a^} encoded in the Weyl function. 



B Appendix: Combinatorial identities 

In the proof of Lcmma l4.1()l we used some combinatorial identities for symmetric 
polynomials. 

Lemma B.l. With the notation of Lemma \4-lU\ as well as that of Defini- 
tion \2.lBl in force, the following identities hold. 

1. With the sums running over all partitions of {1,2, ... ,2k} into disjoint 
subsets I and J with \I\ = \ J\ = k, 

£ AjA^TjTjxj = i £ A?A 2 7 IY, 7 . (B.l) 

i, j i, j 



2. With the sums running over all partitions of {0, 1,2, ... , 2k} into disjoint 
subsets I and J with \I\ — k and \J\ = k + 1, 

AjA'jTjTjX! =± £ A?A*IV, (B.2) 

i,j i,j 

J2 A?A 2 7 r T^) 2 =^ £ A 2 A 2 xjTjj. (B.3) 
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Proof. We give a detailed proof only for the first identity. Almost identical 
technique can be used to verify the other identitities. We start by noting that 
both sides of (|fj.lfl are symmetric polynomials in the variables xi,X2, . . . , X2k- 
In each single variable x a both sides are of degree 3k — 2. Let us denote the 
difference of the two sides by P. We claim that 

1. P is divisible by A(xi, x%, X3, ■ ■ ■ , X2fc)- 

2. P is divisible by A 2 (xi, x 2 , X3, . . . , X2fe). 

Granted Claim 2 we easily conclude that P = and hence the identity is 
proven. Indeed, the degree of P in, say, the x\ variable is 3k — 2, while that of 
A 2 (xi, x 2 , x 3 , x 2 k) is 4fc - 2. 

Furthermore, Claim 2 follows from Claim 1 due to the symmetry of P, the 
skew-symmetry of A(xi, X2, X3, ■ ■ ■ , X2k), arid the fact that any polynomial skew- 
symmetric with respect to all transpositions of variables (in particular, then, the 
polynomial -P/A) must be divisible by A(xi, X2, X3, . . . , x 2 fc). 

Thus it suffices to prove Claim 1. Moreover, by symmetry, it suffices to 
show that P is divisible by x\ — x 2 . The proof proceeds by induction. The 
base case k = is just 1 = 1. Assume the identity holds for k = m — 1. We 
will be evaluating both sides of (|B.1|) at x% = X2. In view of this, the only 
nonvanishing contributions to the left-hand side are coming from sets I and 
J such that I = {1} U V and J = {2} U J', or such that I = {2} U I" and 
J = {1} U J". Thus the only terms contributing to the left-hand side will be 

^2 A {i},7' A /' A {2},j'Aj,r {1}i7 /r//r {2}: j/rj'xix/' 

I', J' 

+ A {2},/"A?//A{ 1}i , / ,/Aj„r {2 } i 7»r/»r {1}i j//r^x2x/». 
i", j" 

Since I' U J' = I" U J" = {3, 4, ... , 2k} =: A, we obtain after evaluating at 
X\ = X2 that the above sums contribute 

2x 2 Af 2M r {2M A ?' A J'I> 
r,j> 

Likewise, following the same steps, the right-hand side contributes 

^ A {2} : A r {2},A Aj,A 2 j,T r j,. 

i', J' 

Since |7'| = | J'| = k — 1, the induction hypothesis ensures that the two con- 
tributions are equal. Thus P vanishes whenever x\ — x 2 , so it is divisible by 
xi — x 2 . This proves Claim 1. □ 

Remark B.2. The only similar identity for symmetric polynomials that we 
have come across in the literature is one involving A(x) 2 , but not F(x), in a 
paper by Anderson on the Toda lattice. 
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